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ABSTRACT 

Age determination is undertaken for nearby early-type (BAF) stars, which constitute attractive 
targets for high-contrast debris disk and planet imaging surveys. Our analysis sequence consists of: 
acquisition of uvbyf3 photometry from catalogs, correction for the effects of extinction, interpolation of 
the photometry onto model atmosphere grids from which atmospheric parameters are determined, and 
finally, comparison to the theoretical isochrones from pre-main sequence through post-main sequence 
stellar evolution models, accounting for the effects of stellar rotation. We calibrate and validate our 
methods at the atmospheric parameter stage by comparing our results to fundamentally determined 
Tgff and logg values. We validate and test our methods at the evolutionary model stage by comparing 
our results on ages to the accepted ages of several benchmark open clusters (IC 2602, a Persei, Pleiades, 
Hyades). Finally, we apply our methods to estimate stellar ages for 3493 field stars, including several 
with directly imaged exoplanet candidates. 

Subject headings: stars: early-type —evolution —fundamental parameters —Hertzsprung-Russell and 
C-M diagrams —planetary systems —astronomical databases: catalogs 


1. INTRODUCTION 

In contrast to other fundamental stellar parameters 
such as mass, radius, and angular momentum - that for 
certain well-studied stars and stellar systems can be an¬ 
chored firmly in observables and simple physics - stellar 
ages for stars other than the Sun have no firm basis. 
Ages are critical, however, for many investigations in¬ 
volving time scales including formation and evolution of 
planetary systems, evolution of debris disks, and inter¬ 
pretation of low mass stars, brown dwarfs, and so-called 
planetary mass objects that are now being detected rou¬ 
tinely as faint point sources near bright stars in high 
contrast imaging surveys. 


1.1. The Era of Direct Imaging of Exoplanets 

Intermediate-mass stars (1.5 — 3.0 Mq) have proven 
themselves attractive targets for planet search work. 
Hints of their importance first arose during initial data 
return from IRAS in the early 1980s, when several A- 
type stars (notably Vega but also (3 Pic and Fomalhaut) 
as well K-star Eps Eri - collectively known as “the fab 
four” - distinguished themselves by showing mid-infrared 
excess emission due to optically thin dust in Kuiper- 
Belt-like locations. Debris disks are signposts of planets, 
which dynamically stir small bodies resulting in dust pro¬ 
duction. Spitzer results in the late 2000s solidified the 
spectral type dependence of debri s disk presence (e.g. 

2006 Wyatt|[2008 ) for stars of common 


Carpenter et al. 


age. For a random sample of field stars, however, the 
primary va riable determining the likelihood of debris is 


stellar age ( [Kains et al. 2011). 

The correlation in radial v elocity studies of giant planet 
frequency with stellar mass (Fischer fc Valenti|2005 Gai- 


dos et al. 2013) is another line ot evidence connecting 


^ NSF Graduate Research Fellow 


planet formation efficiency to stellar mass. The claim is 
that while ^14% of A stars have one or more > lM.T„nit.sr 


companion s at <5 AU, only ^2% of M stars do ( J ohnso n 
et al.||2010 c.f. Lloyd||2M3 Schlaufman &: Winq|2013[ )T 


Consistently interpreted as indicators ot hidden plan- 
ets, debris disks finally had their long-awaited observa¬ 
tional connection to planets with the watershed discovery 
of directly imaged planetary mass companions. These 
were - like the debris disks before them - found first 
around intermediate-mass A-type stars, rather than the 
solar-mass FGK-type stars that had been the subject 
of much observati onal work at high contras t during the 
2000s. HR 8799 (Marois et al. 20081 2010 ) followed by 
Fomalhaut (Kalas et al.|2008() and [3 Pic ( Lagrange et al.| 
2009 2010) have had their planets and indeed one plan¬ 
etary system, digitally captured by ground-based and/or 
space-based high contrast imaging techniques. Of the 
known bona fide planetary mass (< lOMjup) compan¬ 
ions that have been directly imaged, six of the nine are 
located around the three A-type host stars mentioned 
above, with the others associated with lower mass stars 
includ ing the even younger 5-10 Myr old star IRXS 1609- 


2105 (Lafreniere et al. 


2008 


Ireland et al. 2011) and 


brown dwarf 2iVIASS 1207-3933/J[Chauvin^iaL 


2004) 


and the probably older GJ 504 (Kuzuhara et al. 2013). 
Note that to date these directly imaged objects are all 
“super-giant planets” and not solar system giant planet 
analogs (e.g. Jupiter mass or below). 

Based on the early results, the major direct imag¬ 
ing planet searches have attempted to optimize success 
by preferentially observing intermediate-mass, early-type 
stars. The highest masses are avoided due to the lim¬ 
its of contrast. Recent campaigns include those with 
all the major large aperture telescopes: Keck/N1RC2, 
VLT/NACO, Gemini/NlCl, and Subaru/HiCAO. Cur¬ 
rent and near-future campaigns include Project 1640 
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(P1640; Hinkley et al. 2011) at Palomar Observatory, 
Gemini Planet Imager (GPI), operating on the Gemini 
South telescope, VLT/SPHERE, and Subaru/GHARIS. 
The next-generation TMT and E-ELT telescopes both 
fe ature high con t rast i nstruments. 


Mawet et al. (2012) compares instrumental contrast 
curves in their Figure 1. Despite the technological devel¬ 
opments over the past decade, given the as-built contrast 
realities, only the largest, hottest, brightest, and there¬ 
fore the youngest planets, i.e. those less than a few to a 
few hundred Myr in age, are still self-luminous enough to 
be amenable to direct imaging detection. Moving from 
the 3-10 Mjupiter detections at several tens of AU that 
are possible today/soon, to detection of lower mass, more 
Earth-like planets located at smaller, more terrestrial 
zone, separations, will require pushing to higher contrast 
from future space-based platforms. The targets of future 
surveys, whether ground or space, are however not likely 
to be substantially different from the samples targeted in 
today’s ground-based surveys. 

The most important parameter really is age, since the 
brightness of planets decreases so sharply with increasing 
age due to the rapid gravitational contract ion and cooling 
(Eortney et al. 2008 Burrows et al.|2004 ). There is thus 
a premium on identifying the closest, youngest stars. 

1.2. The Age Challenge 

Unlike the other fundamental parameters of stellar 
mass (unambiguously determined from measurements of 
double-line eclipsing binaries and application of Kepler’s 
laws) and stellar radius (unambiguously measured from 
interferometric measurements of angular diameters and 
parallax measurements of distances), there are no di¬ 
rectly interpretable observations leading to stellar age. 

Solar-type stars (~ 0.7—1.4M0, spectral types E6-K5) 
were the early targets of radial velocity planet searches 
and later debris disk searches that can imply the presence 
of planets. Eor these objects, although more work re¬ 
mains to be done, there are established activity-rotation- 
age diagnostics that are driven by the presence of con¬ 
vective outer layers and can serve as p roxies for stellar 
age (e.g. Mamajek fc Hillenbrand|[2008 ). 


Eor stars signihcantiy ditterent from our Sun, how¬ 
ever, and in particular the intermediate-mass stars (^ 
1.5 — 3 .OM 0 , spectral types A0-F5 near the main se¬ 
quence) of interest here, empirical age-dating techniques 
have not been sufficiently established or calibrated. Ages 
have been investigated recently for specific samples of 


Nielsen et al. (|2013): Vigan et al. 

(20121): 

Moor et al. 

imS); SlTit aT-p)0B 

); 

Rhee et af. 

15007 

1 ; 

Lowrance 


Perhaps the most robust ages for young BAE stars 
come from clusters and moving groups, which contain 
not only the early-type stars of interest, but also lower 
mass stars to which the techniques mentioned above can 
be applied. These groups are typically dated using a 
combination of stellar kinematics, lithium abundances, 
rotation-activity indicators, and placement along theo¬ 
retical isochrones in a color-magnitude diagram. The 
statistics of these coeval stellar populations greatly re- 
dnce the uncertainty in derived ages. However, only four 
such groups exist within ~ 60 pc of the Sun and the 
number of early-type members is small. 


Field BAE stars having late-type companions at wide 
separation could have ages estimated using the methods 
valid for F6-K5 age dating. However, these systems are 
not only rare in the solar neighborhood, but consider- 
able effort is require d in establishing companionship e.g . 


Stauffer et al. 


Song et al. HzdOU). Attempts to derive fractional main 


( 1995[);|Barrado y Navascues et al.| (|l997[); 

u 


sequence ages for A-stars bas ed on the evolution of rota- 

ir 


tional velocities are ongoing (Zorec & Royer |2012), but 
this method is undeveloped and a bimodal distribution 
in V sin i for early-type A-stars may inhibit its utility. 
Another method, asteroseismology, which detects low- 
order oscillations in stellar interiors to determine the cen¬ 
tral density and hence age, is a heavily model-dependent 
method, observationally expensive, and best suited for 
older stars with denser cores. 

The most general and quantitative way to age-date 
A0-F5 field stars is through isochrone placement. As 
intermediate-mass stars evolve quickly along the H-R di¬ 
agram, they are better suited for age-dating via isochrone 
placement relative to their low-mass counterparts which 
remai n nearly stationar y on the main sequence for many 
Gyr (Soderblom 2010). Indeed, the mere presence of 


an eariy-type star on the main sequence suggests mod¬ 
erate youth, since the hydrogen burning phase is rela¬ 
tively short-lived. However, isochronal ages are obviously 
model-dependent and they do require precise placement 
of the stars on an H-R diagram implying a parallax. The 


ing metallicitv 1 

Nielsen et al.l 

2013 

), rotation (Collins & 

yinith|iy85 

) and multiplicity (De 

iosa et af.||2014|). 


1.3. Our Approach 

Despite that many nearby BAE stars are well-studied, 
historically, there is no modern data set leading to a set 
of consistently derived stellar ages for this population of 
stars. Here we apply Stromgren photometric techniques, 
and by combining modern stellar atmospheres and mod¬ 
ern stellar evolutionary codes, we develop the methods 
for robust age determination for stars more massive than 
the Sun. The technique uses specific filters, careful cali¬ 
bration, definition of photometric indices, correction for 
any reddening, interpolation from index plots of physi¬ 
cal atmospheric parameters, correction for rotation, and 
finally Bayesian estimation of stellar ages from evolution¬ 
ary models that predict the atmospheric parameters as 
a function of mass and age. 

Specihcally, our work uses high-precision archival 
uvbyP photometry and model atmospheres so as to de¬ 
termine the fundamental stellar atmospheric parameters 
Teff and log 5 . Placing stars accurately in an logTes 
vs. log 5 diagram lea ds to derivation of the ir ages and 
masses. We consider Bressan et al. (2012) evolution¬ 
ary models that include pre-mam sequence evolutionary 
times (2 Myr at 3 Mq and 17 Myr at 1.5 Mq), which are 
a signihcant fraction of any intermediate m ass star’s ab¬ 
solute age, as well as Ekstrom et al. (2012) evolutionary 
models that self-consistentfy account tor stellar rotation, 
which has non-negligible effects on the inferred stellar 
parameters of rapidly rotating early-type stars. Figure 
[2 shows model predictions for the evolution of both phys¬ 
ical and observational parameters. 

The primary sample to which our technique is applied 
in this work consists of 3499 BAE field stars within 100 
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Fig. 1.— Top yane ls: Evolution of logTeff and logg with age for intermediate-mass stars, as predicted by PARSEC evolutionary models 
(|Bressan et al.|2012|l. Bottom panels: Same evolutionary trends for B — V (close to b — y) and My mag, as might be used to discern ages 
from color-magnitude diagram evolution (e.g.|Nielsen et al.|2013|l. While the color and temperature trends reflect one another, the absolute 
magnitude trends are not as strong as the surface gravity trends when the stars are evolving from the main sequence after a few hundred 
Myr. The PARSEC models predict the precision in logp needed to distinguish a 1.5Mq star and a 2.0Mq star evolves from 0.0397 dex at 
~ 30 Myr to 0.0242 dex at 100 Myr to 0.0378 dex at ~ 300 Myr. The precision in log 51 needed to distinguish a I.SMq star and a 3 .OM 0 
evolves from 0.0085 dex at ~ 30 Myr to 0.0694 dex at 100 Myr to 0.5159 dex at ~ 300 Myr. The precision in log^ needed to distinguish a 
2 .OM 0 star and a 3 .OM 0 evolves from 0.0312 dex at ~ 30 Myr to 0.0936 dex at 100 Myr to 0.4781 dex at 300 Myr. 


pc and with uvhy0 photometry available in the Hauck| 
& Mermilliod (19981 catalog, hereafter HM98. The ro¬ 
bustness ot our method is tested at different stages with 
several control samples. To assess the uncertainties in 
our atmospheric parameters we consider (1) 69 stan¬ 
dard s tars f rom iBoyajian et al. (2013) or Napiwotzki 
et al.| (|1993|); (2) 39 double-lined e c lipsin g binaries with 
standard l ogt? from [Torres et al.| (|2010|); (3) 16 other 
stars from Napiwotzld et al.|(|l993), also for examining 
log 5 - To examine isochrone systematics, stars in four 
open clusters are studied (31 members of 1C 2602, 51 
members of a Per, 47 members of the Pleiades, and 47 
members of the Hyades). Some stars belonging to sample 
(1) above are also contained in the large primary sample 
of field stars. 


2. THE STROMGREN PHOTOMETRIC SYSTEM 

Historical use of Strbmgren photometry methods in¬ 
deed has been for the purpose of determining stellar pa¬ 
rameters for e arly-type star s. Recent applicat i ons in ¬ 
clude work by Nieva (|2013|); Dalle Mese et al. (201^); 
Onehag et al. (2009); hlende Prieto & Lambert! (1999). 
An advantage over more traditional color-magnitude di¬ 
agram techniques (Nielsen et al. 2013 De Rosa et al. 
2014) is that distaiice knowledge is not required, so the 
distance-age degeneracy is removed. Also, metallicity ef¬ 
fects are relatively minor (as addressed in an Appendix) 
and rotation effects are well-modelled and can be cor¬ 


rected for (§ 3.3). 


2.1. Description of the Photometric System 



Fig. 2.— The u,v,b,y, H/^^ide H/3narrow passbands. Over¬ 
plotted on an arbi trary scale is the synt hetic spectrum of an AOV 
star generated by |Munari et al.| (|2005|l from an ATLAS9 mode l 
atmosphere. The uvby niter prohles are those of |Besse' m poiTii, 
while the H /3 filter profiles are those originally describe d in 
fordi ||1966|| an d the throughput curves are taken from ICJastelli & 
Kurucz] | |ZU06[ |. 


The uvbyP photometric system is comprised of 
four intermediate -band filters {uvby) first advanced by 
Strbmgren (1 966) plus t he H/ 3 narrow anc^wide filters 
developed by |Crawford| (|1958 ); see Figure [2l Together, 
the two filter sets form a well-calibrated system that was 
specifically designed for studying earlier-type BAF stars, 
for which the hydrogen line strengths and continuum 
slopes in the Balmer region rapidly change with tem¬ 
perature and gravity. 

From the fluxes contained in the six passbands, five 
uvby(3 indices are defined. The color indices, (6 — y) and 
{u — b), and the /3-index, 


f3 = m 


narrow 


H/3viride 5 


( 1 ) 
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are all sensitive to temperature and weakly dependent 
on surface gravity for late A- and F-type stars. The 
Balmer discontinuity index, 

Cl = (w-u) - (n - 5), (2) 

is sensitive to temperature for early type (OB) stars 
and surface gravity for intermediate (AF) spectral types. 
Finally, the metal line index. 


mi = {v-b)- (b-y), 


( 3 ) 


is sensitive to the metallicity [M/H], 

For each index, there is a corresponding intrinsic, 
dereddened index denoted by a naught subscript with e.g 
Co, (6— j/)o, and (u — 6)o, referring to the intrinsic, dered¬ 
dened equivalents of the indices ci, (b—y), and (u — b), re¬ 
spectively. Furthermore, although reddening is expected 
to be negligible for the nearby sources of primary in¬ 
terest to us, automated classification schemes that di¬ 
vide a large sample of stars for analysis into groups cor¬ 
responding to earlier than, around, and later than the 
Balmer maximum will someti mes rely on the reddening- 
indepe ndent indices defined by | Crawford fc Mandwewala| 
(19761 for A-type dwarfs: 


2.88, 0.05 < {b- y)o < 0.22, and (7) F1-G2, IIFV, 
2.60 < P < 2.72, 0.22 < {b—y)o < 0.39. The class values 
in this work were assigned to individual stars based on 
their known spectral types (prov ided in the XHIP cat¬ 
alog (Anderson & Francis 2011)), and /3 values where 
needed. In some instances, Ao- A3 stars assigned to class 
(5) with values of /3 < 2.87, the dereddening procedure 
was unable to proceed. For these cases, stars were either 
assigned to class (1) if they were spectral type AO-Al, or 
to class (6), if they were spectral type A2-A3. 

Depending on the class of an individual star, the 
program then calculates the dereddened indices (b — 
y)o, mo, Co, the color excess E{b—y), brriQ, the absolute V 
magnitude, My, the stellar radius and effective temper¬ 
ature. Notably, the /3 index is unaffected by reddening 
as it is the flux difference between two narrow band fil¬ 
ters with essentially the same central wavelength. Thus, 
no corrections are performed on /3 and this index can be 
used robustly in coarse classification schemes. 

To transform E(b — y) to Ay, we use the extinction 
measurements of|Schlegel et al. (|1998|) and to propagate 
the effects of reddening throug h to the various uvbyp in¬ 
dices w e use the calibrations of ! Crawford &: Mandwewal^ 

(pT^: 


[ci] =ci-0.19(&-y) (4) 

[mi] = mi -I- 0.34(6 — y) (5) 

- 6] = [ci]-I-2[mi]. (6) 


E{mi) =-QME{b - y) (10) 

E{ci) = Q.20E{h - y) (11) 

E{u-b) = lME{h-y). (12) 


Finally, two additional indices useful for early A-type 
stars, ao and r*, are defined as follows: 


ao = 1.36(6 — j/)o -b 0.36mo + 0.18co — 0.2448 (7) 

= (6-j/)o + 0.18[(n-6)o-1.36], (8) 

r* = 0.35ci - 0.07(6 - y) - {P - 2.565). (9) 

Note that r* is a reddening free parameter, and thus 
indifferent to the use of reddened or unreddened photo¬ 
metric indices. 


2.2. Extinction Correction 

Though the sample of nearby stars to which we apply¬ 
ing the Stromgren methodology are assumed to be un- 
extincted or only lightly extincted, interstellar reddening 
is significant for the more distant stars including those 
in the open clusters used in § [^ to test the accuracy of 
the ages derived using our uvEyP methodology. In the 
cases where extinction is thought to be significant, cor¬ 
rections are performed using the UVBYBETl>[jand DEREDE0 
programs for IDL. 

These IDL routines take as input (6 — y), mi, Ci, /3, and 
a class value (between 1-8) that is used to roughly iden¬ 
tify what region of the H-R diagram an individual star 
resides in. For our sample, stars belong to only four 
of the eight possible classes. These classes are summa¬ 
rized as follows: (1) BO-AO, III-V, 2.59 < P < 2.88, 
-0.20 < Co < 1.00, (5) A0-A3, III-V, 2.87 < P < 2.93, 
-0.01 < (6 - y)o < 0.06, (6) A3-F0, III-V, 2.72 < P < 

^ http://idlastro.gsfc.nasa.gov/ftp/pro/astro/uvbybeta. 
pro 

http://idlastro.gsfc.nasa.gov/ftp/pro/astro/deredd.pro 


From these relations, given the intrinsic color index 
(6 — y)o, the reddening free parameters mo, Co, (u — 6)o, 
and ao can be computed. 

In § |4.3| we quantify the effects of extinction and ex¬ 
tinction uncertainty on the final atmospheric parameter 
estimation, Teff,logy. 


2.3. Utility of the Photometric System 

From the four basic Stromgren indices - b — y color, P, 
Cl, and mi - accurate determinations of the stellar atmo¬ 
spheric parameters Teg, logy, and [M/E[] are possible for 


Crawford||19791 ILester et al. 

119861 |01sen||198 

81 ISmallev 

I993| Smalley & Dworetsky 

19951 Clem et al.||2004D, or 

theoretical (e.g. Balona||1984 

1 |Moon & Dworetsky||1985 

Napiwotzki et al.|1993l|Balona|1994l|Leieune et al.| 


Castelli & Kurucz 

2006 

2004| 

Onehag et al. 

2009 

) cal- 


in Teff are claimed as a chiev able and we reassess these 
uncertainties ourselves § |4.3| 


3. DETERMINATION OF ATMOSPHERIC 
PARAMETERS Teff, log G 

3.1. Procedure 

Once equipped with uvbyp colors and indices and un¬ 
derstanding the effects of extinction, arriving at the fun¬ 
damental parameters Teg and logy for program stars, 
proceeds by interpolation among theoretical color grids 
(generated by convolving filter sensitivity curves with 
model atmospheres) or explicit formulae (often polyno¬ 
mials) that can be derived empirically or using the theo¬ 
retical color grids. In both cases, calibration to a sample 
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Fig. 3.— Top: Three relevant uvbyfi spaces for atmospheric parameter determination of our sample of BAF stars with uvbyfi photometry 
in the HM98 catalog, and located within 100 pc of the Sun. Two stars were excluded from these figures for favorable scaling: Castor, which 
is an outlier in all three planes (/3 <2.4, ao > 1, b — v >0.6~). and HD 17300, a poorly studied F3V star with b — y >0.6. Bottom: The same 
plots as above, with the model color grids of|Castelli &: Kurucz| (|2006l |2004|| overlaid in the relevant regions of parameter space. The lines 
of constant Teg (largely vertical) and of constant logg (largely horizontal) are annotated with their corresponding values. Some outliers 
have been pruned, and irrelevant groups of stars eliminated, for clarity in this second plot. 


of stars with atmospheric parameters that have been in¬ 
dependently deter mined through fundam ental physics is 
required. See e.g. Figueras et al. (19911 for further de¬ 
scription. 

Numerous calibrations, both theoretical and empirical, 
of the uvby B photometric syst e m exist. Fo r this work 
we use the Castelli & Kurucz (2006 2004) color grids 
generated from solar metallicity (Z=0.017, in this case) 
ATLAS9 model atmospheres using a microturbulent ve¬ 
locity parameter of ^ = 0 km s“^ and the new ODF. We 
do not use the alpha-enhanced color grids. The CTids are 
readily available from F. Castellj^or R. Kuruc^ 

Prior to assigning atmospheric parameters to our pro¬ 
gram stars directly from the model grids, we first inves¬ 
tigated the accuracy of the models on samples of BAF 
stars with fundamentally determined T^s (through in¬ 
terferometric measurements of the angular diameter and 
estimations of the total integrated flux) and log 5 (from 
measurements of the masses and radii of double lined 
eclipsing bina ries). We describe these validation proce¬ 
dures in § |4.1| and § |4.2[ 

Atmospheric parameter determination occurs in three 
different observational Stromgren planes depending on 
the temperature regime (see Figure |^; this is in order 
to avoid the degeneracies that are present in all single 
observational planes when mapped onto the physical pa¬ 
rameter space of logTeff and logi/. _ 

Building o ff of the original work of e.g. IStromg: 

and later |Na 


([I^[T9^ 


ren 


Moon & Dworetsky (1985) 


'lapi- 


wotzki et ai. (1993), suggested assigning physical param- 


eters in the toliowing three regimes: for cool stars (Te® < 

http: //wwwuser. oats. inaf . it/castelli 
® http://kurucz.harvard.edu/grlds/gridPOOODFNEW/ 

uvbypOOkOodfnew.dat 


8500 K), P or {b—y) can be used as a temperature indica¬ 
tor and Co is a surface gravity indicator; for intermediate 
temperature stars (8500 K < < 11000 K), the tem¬ 

perature indicator is oq and surface gravity indicator r*; 
finally, for hot stars (Teff > 11000 K), the cq or the [u — b] 
indices can be used as a temperature indicator while $ 
is a gravity indicator (note that the role of /3 is reversed 
for hot stars compared to its role for cool stars). We 
adopt here Ci vs. /3 for the hottest stars, Qq vs. r* for 
the intermediate temperatures, and {b— y) vs. ci for the 
cooler stars. 

Choosing the appropriate plane for parameter deter¬ 
mination effectively means establishing a crude temper¬ 
ature sequence prior to fine parameter determination; in 
this, the P index is critical. Because the P index switches 
from being a temperature indicator to a gravity indicator 
in the temperature range of interest to us (spectral type 
B0-F5, luminosity class IV /V stars), atmospheric param¬ 
eter determination proceeds depending on the tempera¬ 
ture regime. For the T^g and log 5 calibrations described 
below, temperature information existed for all of the cali¬ 
bration stars, though this is not the case for our program 
stars. In the general case we must rely on photometric 
classification to assign stars to the late, intermediate, and 
early groups, and then proceed to determine atmospheric 
p aramet e rs in t he relevant uvbyP planes. 


Moon (1985) provides a scheme, present in the 
UVBYBETA lUL routine, for roughly identifying the region 
of the H-R diagram in which a star resides. However, 
because our primary sample of held stars are assumed 
to be unextincted, and because the UVBYBETA program 
relies on user-inputted class values based on unverified 
spectral types from the literature, we opt for a classifica¬ 
tion scheme based solely on the uvbyP photometry. 
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Monguio et al. (2014), hereafter M14, designed a so- 
phisticated classifi cation scheme, based on the work of 
Stromgren (1966). The M14 scheme places stars into 
early (HU-Al)j, intermediate (AO-A3), and late (later 
than A3) groups based solely on /3, the reddened color 
(b—y), and the reddening-free parameters [ci], [m-i], [u — 
b]. The M14 scheme imp roves upon the previous method 
of Figueras et al. (1991) by imposing two new conditions 
(see their h'igure 2 tor the complete scheme) intended 
to prevent the erroneous classification of some stars. For 
our sample of 3499 field stars (see §[^, there are 699 stars 
lacking /3 photometry, all but three of which cannot be 
classified by the M14 scheme. For such cases, we rely on 
supplementary spectral type information and manually 
assign these unclassified stars to the late group. Using 
the M14 scheme, the final makeup of our field star sample 
is 85.9% late, 8.4% intermediate, and 5.7% early. 

3.2. Sample and Numerical Methods 

For all stars i n this work, uvby0 photome try is ac¬ 
quired from the |Hauck & Mermilliod (1998) compila¬ 
tion (hereafter HiV198), unless otherwise noted. HM98 
provides the most extensive compilation of uvby/3 pho¬ 
tometric measurements, taken from the literature and 
complete to the end of 1996 (the photometric system 
has seen less frequent usage/publication in more modern 
times). The HM98 compilation includes 105,873 individ¬ 
ual photometric measurements for 63,313 different stars, 
culled from 533 distinct sources, and are presented both 
as individual measurements and weighted means of the 
literature values. 

The HM98 catalog provides (6 — y),mi,ci, and /3 and 
the associated errors in each parameter if available. From 
these indices oq and r* are computed according to Equa¬ 
tions (7), ( 8 ) & (9). The ATLAS9 uvbyP grids provide 
a means of translating from (6 — y, mi, Ci, /3, Oq, r*) to a 
precise combination of (Teff, log <?). Interpolation within 
the model grids is performed on the appropriate grid: 
{{b — y) vs. Cl for the late group, oq vs. r* for the 
intermediate group, and ci vs. /3 for the early group). 

The interpolation is linear and performed using the 
SciPy routine griddata. Importantly, the model log g 
values are first converted into linear space so that g is 
determined from the linear interpolation procedure be¬ 
fore being brought back into log space. The model grids 
used in this work are spaced by 250 K in Teg and 0.5 dex 
in log 5 . To improve the precision of our method of atmo¬ 
spheric parameter determination in the future, it would 
be favorable to use model color grids that have been cal¬ 
culated at finer resolutions, particularly in log 5 , directly 
from model atmospheres. However, the grid spacings 
stated above are fairly standardized among extant uvbyfi 
grids. 

3.3. Rotational Velocity Correction 

Early-type stars are rapid rotators, with rotational ve¬ 
locities of vsini > 150 km s“^ being typical. For a ro¬ 
tating star, both surface gravity and effective tempera¬ 
ture decrease from the poles to the equator, changing the 
mean gravity and te mperature of a rapid rotator relative 
to a slower rotator (Sweet & Roy 1953). Vega, rotat¬ 


ing with an inferred equatorial velocity of Veq 270 km 
s“^ at a nearly pole-on inclination, has measured pole- 
to-equator gradients in Teff and logg that are ~ 2400 K 



Fig. 4. — Vectors showing the magnitude and direction of the 
rotational velocity corrections at 100 (black), 200, and 300 (light 
grey) km s~^ for a grid of points in log(Teff)-log^ space, with 
PARSEC isochrones overlaid for reference. While typical A-type 
stars rotate at about 150 km s“^, high-contrast imaging targets 
are sometimes selected for slow rotation and hence favorable incli¬ 
nations, typically usinz <50 km s“^ or within the darkest black 
vectors. For rapid rotators, a 100% increase in the inferred age due 
to rotational effects is not uncommon. 


and ^ 0.5 dex, respectively ( Peterson et al.|[2006 ). The 
apparent luminosity change due to rotation depends on 
the inclination: a pole-on (i = 0 °) rapid rotator appears 
more luminous than a nonrotating star of the same mass, 
while an edge-on {i = 90°) rapid rotator appears less lu¬ 


minous than a nonrotating star of the same mass. Sweet 


& Roy (1953) found that a (usini)^ correction factor 


could describe the changes in luminosity, gravity, and 
temperature. 

The net effect of stellar rotation on inferred age is to 
make a rapid rotator appear cooler, more luminous, and 
hence older when compared to a nonrotating star of the 
same mass (or more massive when compared to a non¬ 
rotating star of the same age). Optical colors can be 
affected since the spectral lines of early type stars are 
strong and broad. Kraft & Wrubel (1965) demonstrated 
specifically in the Stromgren system that the effects are 
predominantly in the gravity indicators (ci, which then 
also affects the other gravity indicator r*) and less so in 
t he temperature indicator s {b — y, which then affects oq). 


Figueras & Blasi (1998), hereafter FB98, used Monte- 
Carlo simulations to investigate the effect of rapid ro¬ 
tation on the measured uvbyjd indices, derived atmo¬ 
spheric parameters, and hence isochronal ages of early- 
type stars. Those authors concluded that stellar rotation 
conspires to artificially enhance isochronal ages derived 
through uvby/3 photometric methods by 30-50% on aver¬ 
age. 

To mitigate the effect of stellar rotation on the pa¬ 
rameters Teff and log(( 7 ), FB98 presented the following 
corrective formulae for stars with Teff > 11000 K: 


ATeff = 0.0167(u sin i)^ + 218, (13) 

Alog 5 = 2.10 X 10"®(usini)2 -t 0.034. (14) 

For stars with 8500K < Teff < llOOOK, the analogous 
formulae are: 
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ATeff = 0.0187(z; sin + 150, (15) 

Alogg = 2.92 x 10-®(usini)2 + 0.048. (16) 


In both cases, ATeff and Alogg are added to the Teff 
and logg values derived from uvbyfi photometry. 

Notably, the rotational velocity correction is dependent 
on whether the star belongs to the early, intermediate, 
or late group. Specifically, FB98 define three regimes: 
Teff <8830 K (no correction), 8830 K<Teff<9700 K (cor¬ 
rection for intermediate AO-A3 stars), Teff>9700 K (cor- 
r ection for stars ear lier than A3). 


Song et al. (2001), who performed a similar isochronal 
age analysis of A-type stars using uvbyjd photometry, 
extended the FB98 rotation corrections to stars earlier 
and later than B7 and A4, respectively. In the present 
work, a more conservative approach is taken and the ro¬ 
tation correction is applied only to stars in the early or 
intermediate groups, as determined by the classification 


scheme discussed in § |3.1 
tified by the abundance o: 


This decision was partly jus- 
late-type stars that fall below 


the ZAMS in the open cluster tests (§ct, for which the 
rotation correction would have a small ^ue to the lower 
rotational velocities of late-type stars) but exacerbating 
effect on these stars whose surface gravities are already 
thought to be overestimated. 

We include these corrections and, as illustrated in Fig¬ 
ure 1^ emphasize that in their absence we would err on 
the side of over-estimating the age of a star, meaning 
conservatively overestimating rather than underestimat¬ 
ing companion masses based on assumed ages. As an 
example, for a star with Tes « 13,275 K and logg « 4.1, 
assumed to be rotating edge-on at 300 km s“^, neglect¬ 
ing to apply the rotation correction would result in an 
age of ^ 100 Myr. Applying the rotation correction to 
this star results in an age of ~ 10 Myr. 

Of note, the FB98 corrections were derived for atmo¬ 
spheric para meters determined using the synthetic uvbyfd 
color grids of|Moon & Dworetsky| (|1985|) . It is estimated 
that any differences in derived atmospheric parameters 
resulting from the use of c olor grids other than those of 


Moon & Dworetsky 


surement errors m t 


(1985) are less than t he t ypical mea- 
lose parameters. In § 4.3 we quantify 
the effects of rotation and rotation correction uncertainty 
on the final atmospheric parameter estimation, Teff, log g. 


4. CALIBRATION AND VALIDATION USING THE 
HM98 CATALOG 


In this section we assess the effective temperatures and 
surface gravities derived from atmospheric models and 
uvbyP color gri ds re lative to fundamentally dete rmined 
temperatures (§ 4.1) and surface gravities (§ 4.2). 


4.1. Effective Temperature 

A fundamental determination of Teff is possible 
through an interferometric measurement of the stellar 
angular diameter and an estimate of the total integrated 
flux. We gathered 69 stars (listed in Table with funda¬ 
mental Teff measurements from the literature and deter¬ 
mine photometric temperatures for these objects from 
interpolation of uvbyjd photometry in ATLAS9 model 
grids. 


Fundame ntal Teff values wer e sourced from |Boyajian| 
et al. (2013), hereafter B13, and Napiwotzki et ^ ( [1993^ , 


hereafter N93. Several stars have multiple interferomet- 
ric measurements of the stellar radius, and hence mul¬ 
tiple fundamental Teff determinations. For these stars, 
identified as those objects with multiple radius references 
in Table the mean Teff and standard deviation were 
taken as the fundamental measurement and standard er¬ 
ror. Among the 16 stars with multiple fundamental Teff 
determinations by between 2 and 5 authors, there is a 
scatter of typically several percent (with 0.1-4% range). 

Additional characteristics of the Tgff “standard” stars 
are summarized as follows: spectral types B0-F9, lumi¬ 
nosity classes III-V, 2 km s“^ < usini < 316 km s“^, 
mean and median vsini of 58 and 26 km s“^, respec¬ 
tively, 2.6 pc < d < 493 pc, and a mean and median 
[Fe/H] of -0.08 and -0.06 dex, respectively. Line-of-sight 
rotational velocit ies were acquired from the|Glebocki fc| 


Gnacinski (2005|) compilation and [Fe/H] values were 
taken frorn SliVlBAD. Variability and multiplicity were 
considered, and our sample is believed to be free of any 
possible contamination due to either of these effects. 

From the HM98 compilation we retrieved uvbyfd pho¬ 
tometry for these “effective temperature standards.” 
The effect of reddening was considered for the hot¬ 
ter, statistically more distant stars in the N93 sample. 
Gomparing mean uvhyjd photometry from HM98 with 
the dereddened photometry presented in N93 revealed 
that nearly all of these stars have negligible reddening 
{E{b — y) < 0.001 mag). The exceptions are HD 82328, 
HD 97603, HD 102870, and HD 126660 with color ex¬ 
cesses of E{b— y) = 0.010, 0.003, 0.011, and 0.022 mag, 
respectively. Inspection of Table Vindicates that despite 
the use of the reddened HM98 pnotometry the Tgff de¬ 
terminations for three of these four stars are still of high 
accuracy. For HD 97603, there is a discrepancy of > 300 
K between the fundamental and photometric tempera¬ 
tures. However, the uvbyfi Tgg using reddened photom¬ 
etry for this star is actually hotter than the fundamental 
Tgff. Notably, the author-to-author dispersion in multi¬ 
ple fundamental Tgff determinations for HD 97603 is also 
rather large. As such, the HM98 photometry was deemed 
suitable for all of the “effective temperature standards.” 

For the sake of completeness, differen t model color 
grids were inv estigated, including those of Fitzpatrick fc| 


Massa (2005), which were recently calibrated t or early 
group stars, and those of|6nehag et al. (2009), which 
were calibrated from MARGb model atrnospfieres for 
stars cooler than 7000 K. We found the grids that best 
matched the fundamental effective temperatures were 
the ATLAS9 grids of solar metallicity with no alpha- 
enhancement, microturbulent velocity of 0 km s“^, and 
using the new opacity distribution function (ODF). The 
ATLAS9 grids with microturbulent velocity of 2 km s“^ 
were also tested, but were found to worsen both the frac¬ 
tional Tgff error and scatter, though only nominally (by 
a few tenths of a percent). 

For the early group stars, temperature determinations 
were attempted in both the ci — /3 and [m — b]-P planes. 
The Cl index was found to be a far better temperature 
indicator in this regime, with the [u — 6] index underesti¬ 
mating Tgff relative to the fundamental values >10% on 
average. Temperature determinations in the Ci —/3 plane, 
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Fig. 5. — Top: Comparison of the temperatures derived from the 
ATLAS9 uvbyfi color grids (Tuvby) S'lid the fundamental effective 
temperatures taken from B13 and N93. Bottom: Ratio of 

uvbyp temperature to fundamental temperature, as a function of 
Fuvby For the majority of stars, the uvbyf} grids can predict Tgir 
to within ~ 5% without any additional correction factors. 


however, were only « 1.9% cooler than the fundamental 
values, regardless of whether ci or the dereddened index 
Co was used. This is not surprising as the Ci — /3 plane is 
not particularly susceptible to reddening. 

At intermediate temperatures, the a^—r* plane is used. 
In this regime, the ATLAS9 grids were found to overesti¬ 
mate Teff by « 2.0% relative to the fundamental values. 

Finally, for the late group stars temperature determi¬ 
nations were attempted in the (b — y) — Ci and P — Ci 
planes. In this regime, {b — y) was found to be a su¬ 
perior temperature indicator, improving the mean frac¬ 
tional error marginally and reducing the RMS scatter by 
more than 1%. In this group, the model grids overpre¬ 
dict Teff by « 2.4% on average, regardless of whether the 
reddened or dereddened indices are used. 

Figure shows a comparison of the temperatures de¬ 
rived from the ATLAS9 uvbyP color grids and the funda¬ 
mental effective temperatures given in B13 and N93. For 
the majority of stars the color grids can predict the effec¬ 
tive temperature to within about 5 %. A slight system¬ 
atic trend is noted in Figure such that the model color 
grids overpredict Teff at low temperatures and underpre¬ 
dict Teff at high temperatures. We attempt to correct 
for this systematic effect by applying Teff offsets in three 
regimes according to the mean behavior of each group: 
late and intermediate group stars were shifted to cooler 
temperatures by 2.4% and 2.0%, respectively, and early 
group stars were shifted by 1.9% toward hotter tempera¬ 
tures. After offsets were applied, the remaining RMS er¬ 
ror in temperature determinations for these “standard” 
stars was 3.3%, 2.5%, and 3.5% for the late, intermedi¬ 


ate, and early groups, respectively, or 3.1% overall. 

Taking the uncertainties or dispersions in the funda¬ 
mental Teff determinations as the standard error, there 
is typically a 5-6 a discrepancy between the fundamen¬ 
tal and photometric Tgff determinations. However, given 
the large author-to-author dispersion observed for stars 
with multiple fundamental Teff determinations, it is likely 
that the formal errors on these measurements are un¬ 
derestimated. Notably, N93 does not publish errors for 
the fundamental Teff values, which are literature means. 
However, those authors did find fractional errors in their 
photometric Teff ranging from 2.5-4% for BA stars. 

In §[^ we opted not to apply systematic offsets, instead 
assigning Teff uncertainties in three regimes according to 
the average fractional uncertainties noted in each group. 
In our final Teff determinations for our field star sample 
(§ [T]) we attempted to correct for the slight temperature 
sy^ematics and applied offsets, using the magnitude of 
the remaining RMS error (for all groups considered col¬ 
lectively) as the domin ant s ource of uncertainty in our 
Teff measurement (see § |4.3[ ). 

As demonstrated in Figure rotational effects on our 
temperature determinations for the Tgff standards were 
investigated. Notably, the FB98 usini corrections ap¬ 
pear to enhance the discrepancy between our tempera¬ 
ture determinations and the fundamental temperatures 
for the late and intermediate groups, while moderately 
improving the accuracy for the early group. For the late 
group this is expected, as the correction formulae were 
originally derived for intermediate and early group stars. 
Notably, however, only two stars in the calibration sam¬ 
ple exhibit projected rotational velocities > 200 km s“^. 
We ex amine the utility of the v sin i correction further in 

§lll&§i 


The efiect of metallicity on the determination of Tgff 
from the uvbyP grids is investigated in Figure show¬ 
ing the ratio of the grid-determined temperature to the 
fundamental temperature as a function of [Fe/H]. The 
sample of temperature standards spans a large range in 
metallicity, yet there is no indication of any systematic 
effect with [Fe/H], justifying our choice to assume solar 
metallicity throughout this work (see further discussion 
of metallicity effects in the Appendix). 

The effect of reddening on our temperature deter¬ 
minations was considered but since the vast majority 
of sources with fundamental effective temperatures are 
nearby, no significant reddening was expected. Indeed, 
no indication of a systematic trend of the temperature 
residuals as a function of distance was noted. 

In summary our findings that the ATLAS9 predicted 
Teff values are ^ 2% hotter than fundament al values 
for AF star s are consistent with the results of IBertonel 
(20041, who found 4-8% shifts warmer in Teff from 
ATLAS9 models to spectrophotometry relative to 
Teff values determined from the infrared flux method 
(IRFM). We attempt systematic corrections with offsets 
of magnitude ~ 2% according to group, and the remain¬ 
ing RMS error between uvbyP temperatures and funda¬ 
mental values is ^ 3%. 


et al. 
tits o: 


4.2. Surface Gravity 

To assess the surface gravities derived from the uvbyP 
grids, we compare to results on both double-lined eclips¬ 
ing binary and spectroscopic samples. 
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TABLE 1 

Stars with fundamental determinations of Tgff through Interferometry 


HD 

Sp. Type 

Tfund 

(K) 

Radius Ref.“ 

'^uvbylS 

(K) 

Quvby^ 

(dex) 

[Fe/H] 

(dex) 

V sin i 
(km s~^) 

(b-y) 

(mag) 

mi 

(mag) 

Cl 

(mag) 

/3 

(mag) 

4614 

F9V 

5973 ± 8 

3 

5915 

4.442 

-0.28 

1.8 

0.372 

0.185 

0.275 

2.588 

5015 

F8V 

5965 ± 35 

3 

6057 

3.699 

0.04 

8.6 

0.349 

0.174 

0.423 

2.613 

5448 

A5V 

8070 

18 

8350 

3.964 

— 

69.3 

0.068 

0.189 

1.058 

2.866 

6210 

F6Vb 

6089 ± 35 

1 

5992 

3.343 

-0.01 

40.9 

0.356 

0.183 

0.475 

2.615 

9826 

F8V 

6102 ± 75 

2,4 

6084 

3.786 

0.08 

8.7 

0.346 

0.176 

0.415 

2.629 

16765 

F7V 

6356 ± 46 

1 

6330 

4.408 

-0.15 

30.5 

0.318 

0.160 

0.355 

2.647 

16895 

F7V 

6153 ± 25 

3 

6251 

4.118 

0.00 

8.6 

0.325 

0.160 

0.392 

2.625 

17081 

B7V 

12820 

18 

12979 

3.749 

0.24 

23.3 

-0.057 

0.104 

0.605 

2.717 

19994 

F8.5V 

5916 ± 98 

2 

5971 

3.529 

0.17 

7.2 

0.361 

0.185 

0.422 

2.631 

22484 

F9IV-V 

5998 ± 39 

3 

5954 

3.807 

-0.09 

3.7 

0.367 

0.173 

0.376 

2.615 

30652 

F6IV-V 

6570 ± 131 

3,6 

6482 

4.308 

0.00 

15.5 

0.298 

0.163 

0.415 

2.652 

32630 

B3V 

17580 

18 

16536 

4.068 

— 

98.2 

-0.085 

0.104 

0.318 

2.684 

34816 

B0.5IV 

27580 

18 

28045 

4.286 

-0.06 

29.5 

-0.119 

0.073 

-0.061 

2.602 

35468 

B2III 

21230 

18 

21122 

3.724 

-0.07 

53.8 

-0.103 

0.076 

0.109 

2.613 

38899 

B9IV 

10790 

18 

11027 

3.978 

-0.16 

25.9 

-0.032 

0.141 

0.906 

2.825 

47105 

AOIV 

9240 

18 

9226 

3.537 

-0.28 

13.3 

0.007 

0.149 

1.186 

2.865 

48737 

F5IV-V 

6478 ±21 

3 

6510 

3.784 

0.14 

61.8 

0.287 

0.169 

0.549 

2.669 

48915 

AOmAlVa 

9755 ± 47 

7,8,9,10,11 

9971 

4.316 

0.36 

15.8 

-0.005 

0.162 

0.980 

2.907 

49933 

F2Vb 

6635 ± 90 

12 

6714 

4.378 

-0.39 

9.9 

0.270 

0.127 

0.460 

2.662 

56537 

A3Vb 

7932 ± 62 

3 

8725 

4.000 

— 

152 

0.047 

0.198 

1.054 

2.875 

58946 

FOVb 

6954 ± 216 

3,18 

7168 

4.319 

-0.25 

52.3 

0.215 

0.155 

0.615 

2.713 

61421 

F5IV-V 

6563 ± 24 

11,13,14,15,18 

6651 

3.983 

-0.02 

4.7 

0.272 

0.167 

0.532 

2.671 

63922 

BOIII 

29980 

18 

29973 

4.252 

0.16 

40.7 

-0.122 

0.043 

-0.092 

2.590 

69897 

F6V 

6130 ± 58 

1 

6339 

4.290 

-0.26 

4.3 

0.315 

0.149 

0.384 

2.635 

76644 

A7IV 

7840 

18 

8232 

4.428 

-0.03 

142 

0.104 

0.216 

0.856 

2.843 

80007 

A2IV 

9240 

18 

9139 

3.240 

— 

126 

0.004 

0.140 

1.273 

2.836 

81937 

FOIVb 

6651 ± 27 

3 

7102 

3.840 

0.17 

146 

0.211 

0.180 

0.752 

2.733 

82328 

F5.5IV-V 

6299 ± 61 

3,18 

6322 

3.873 

-0.16 

7.1 

0.314 

0.153 

0.463 

2.646 

90839 

F8V 

6203 ± 56 

3 

6145 

4.330 

-0.11 

8.6 

0.341 

0.171 

0.333 

2.618 

90994 

BOV 

14010 

18 

14282 

4.219 

— 

84.5 

-0.066 

0.111 

0.466 

2.730 

95418 

AlIV 

9181 ± 11 

3,18 

9695 

3.899 

-0.03 

40.8 

-0.006 

0.158 

1.088 

2.880 

97603 

A5IV(n) 

8086 ± 169 

3,6,18 

8423 

4.000 

-0.18 

177 

0.067 

0.195 

1.037 

2.869 

102647 

A3Va 

8625 ± 175 

5,6,18 

8775 

4.188 

0.07 

118 

0.043 

0.211 

0.973 

2.899 

102870 

F8.5IV-V 

6047 ± 7 

3,18 

6026 

3.689 

0.12 

5.4 

0.354 

0.187 

0.416 

2.628 

118098 

A2Van 

8097 ± 43 

3 

8518 

4.163 

-0.26 

200 

0.065 

0.183 

1.006 

2.875 

118716 

BlIII 

25740 

18 

23262 

3.886 

— 

113 

-0.112 

0.058 

0.040 

2.608 

120136 

F7IV-V 

6620 ± 67 

2 

6293 

3.933 

0.24 

14.8 

0.318 

0.177 

0.439 

2.656 

122408 

A3V 

8420 

18 

8326 

3.500 

-0.27 

168 

0.062 

0.164 

1.177 

2.843 

126660 

F7V 

6202 ± 35 

3,6,18 

6171 

3.881 

-0.02 

27.7 

0.334 

0.156 

0.418 

2.644 

128167 

F4VkF2mFl 

6687 ± 252 

3,18 

6860 

4.439 

-0.32 

9.3 

0.254 

0.134 

0.480 

2.679 

130948 

F9IV-V 

5787 ± 57 

1 

5899 

4.065 

-0.05 

6.3 

0.374 

0.191 

0.321 

2.625 

136202 

F8IV 

5661 ± 87 

1 

6062 

3.683 

-0.04 

4.9 

0.348 

0.170 

0.427 

2.620 

141795 

kA2hA5mA7V 

7928 ± 88 

3 

8584 

4.346 

0.38 

33.1 

0.066 

0.224 

0.950 

2.885 

142860 

F6V 

6295 ± 74 

3,6 

6295 

4.130 

-0.17 

9.9 

0.319 

0.150 

0.401 

2.633 

144470 

BIV 

25710 

18 

25249 

4.352 

— 

107 

-0.112 

0.043 

-0.005 

2.621 

162003 

F5IV-V 

5928 ± 81 

3 

6469 

3.916 

-0.03 

11.9 

0.294 

0.147 

0.497 

2.661 

164259 

F2V 

6454 ± 113 

3 

6820 

4.121 

-0.03 

66.4 

0.253 

0.153 

0.560 

2.690 

168151 

F5Vb 

6221 ± 39 

1 

6600 

4.203 

-0.28 

9.7 

0.281 

0.143 

0.472 

2.653 

169022 

B9.5III 

9420 

18 

9354 

3.117 

— 

196 

0.016 

0.102 

1.176 

2.778 

172167 

AOVa 

9600 

18 

9507 

3.977 

-0.56 

22.8 

0.003 

0.157 

1.088 

2.903 

173667 

F5.5IV-V 

6333 ± 37 

3,18 

6308 

3.777 

-0.03 

16.3 

0.314 

0.150 

0.484 

2.652 

177724 

AOIV-Vnn 

9078 ± 86 

3 

9391 

3.870 

-0.52 

316 

0.013 

0.146 

1.080 

2.875 

181420 

F2V 

6283 ± 106 

16 

6607 

4.187 

-0.03 

17.1 

0.280 

0.157 

0.477 

2.657 

185395 

F3+V 

6516 ± 203 

3,4 

6778 

4.296 

0.02 

5.8 

0.261 

0.157 

0.502 

2.688 

187637 

F5V 

6155 ± 85 

16 

6192 

4.103 

-0.09 

5.4 

0.333 

0.151 

0.380 

2.631 

190993 

B3V 

17400 

18 

16894 

4.195 

-0.14 

140 

-0.083 

0.100 

0.295 

2.686 

193432 

B9.5V 

9950 

18 

10411 

3.928 

-0.15 

23.4 

-0.021 

0.134 

1.015 

2.852 

193924 

B2IV 

17590 

18 

17469 

3.928 

— 

15.5 

-0.092 

0.087 

0.271 

2.662 

196867 

B9IV 

10960 

18 

10837 

3.861 

-0.06 

144 

-0.019 

0.125 

0.889 

2.796 

209952 

B7IV 

13850 

18 

13238 

3.913 

— 

215 

-0.061 

0.105 

0.576 

2.728 

210027 

F5V 

6324 ± 139 

6 

6496 

4.187 

-0.13 

8.6 

0.294 

0.161 

0.446 

2.664 

210418 

A2Vb 

7872 ± 82 

3 

8596 

3.966 

-0.38 

136 

0.047 

0.161 

1.091 

2.886 

213558 

AlVb 

9050 ± 157 

3 

9614 

4.175 

— 

128 

0.002 

0.170 

1.032 

2.908 

215648 

F6V 

6090 ± 22 

3 

6198 

3.950 

-0.26 

7.7 

0.331 

0.147 

0.407 

2.626 

216956 

A4V 

8564 ± 105 

5,18 

8857 

4.198 

0.20 

85.1 

0.037 

0.206 

0.990 

2.906 

218396 

F0+(A Boo) 

7163 ± 84 

17 

7540 

4.435 

— 

47.2 

0.178 

0.146 

0.678 

2.739 

219623 

F8V 

6285 ± 94 

1 

6061 

3.85 

0.04 

4.9 

0.351 

0.169 

0.395 

2.624 

222368 

F7V 

6192 ± 26 

3 

6207 

3.988 

-0.14 

6.1 

0.330 

0.163 

0.399 

2.625 

222603 

A7V 

7734 ± 80 

1 

8167 

4.318 

— 

62.8 

0.105 

0.203 

0.891 

2.826 


Interferometric radii references: ( 1 ) Bovaiian et al.lll 201 3 ll 

— HfegiTf 


I Di Fol co et all || 2004 ||, ( 6 )lvan tse lie ^ ^on^raunllf 
(lU) [Kervella et al.lJ^UU^j 


|et al. 
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Kervella et al. 
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Fig. 6 .— Ratio of the uvby^ temperature to fundamental tem¬ 
perature as a function of usini, for the late (left), intermediate 
(middle), and early (right), group stars. The solid horizontal 
colored lines indicate the mean ratios in each case. The arrows 
reperesent both the magnitude and direction of change to the ra¬ 
tio Tuy}jy /T^y^Y^d after applying the FB98 rotation corrections. The 
dashed horizontal colored lines indicate the mean ratios after ap¬ 
plication of the rotation correction. The rotation correction ap¬ 
pears to improve temperature estimates for early group stars, but 
worsen estimates for the late and intermediate groups. Notably, 
however, the vast majority of Teff standards are slowly rotating 
(usinz < 150 km s“^). Note one rapidly rotating intermediate 
group star extends beyond the scale of the figure, with a rotation 
corrected T^vby /'^fund ratio of Pi 1.26. 


• 

• 

M * 

-r >1 


late 

intermediate 

early 


- 0.6 - 0.4 - 0.2 0.0 0.2 0.4 

[Fe/HJ 


Fig. 7.— Ratio of the uvby^ temperature to fundamental tem¬ 
perature as a function of [Fe/H]. There is no indication that the 
grids systematically overestimate or underestimate Tgff for differ¬ 
ent values of [Fe/H]. 


4.2.1. Comparison with Double-Line Eclipsing Binaries 


Torres et al.| (2010|l compiled an extensive catalog of 
95 double-lined eclipsing binaries with fundamentally de¬ 
termined surface gravities for all 180 individual stars. 
Eclipsing binary systems allow for dynamical determi¬ 
nations of the component masses and geometrical deter¬ 
minations of the component radii. From the mass and 
radius of an individual component, the Newtonian sur¬ 
face gravity, g = GM/E? can be calculated. 

From these systems, 39 of the primary components 
have uvhyl3 photometry available for determining sur¬ 
face gravities using our methodology. The spectral type 
range for these systems is 08-F2, with luminosity classes 
of IV, V. The mass ratio (primary/secondary) for these 
systems ranges from « 1.00-1.79, and the orbital peri¬ 
ods of the primaries range from « 1.57-8.44 days. In the 
cases of low mass ratios, the primary and secondary com¬ 
ponents should have nearly identical fundamental param¬ 
eters, assuming they are coeval. In the cases of high mass 
ratios, given that the individual components are presum¬ 
ably unresolved, we assume that the primary dominates 
the uvby(j photometry. For both cases (of low and high 
mass ratios), we assume that the photometry allows for 
accurate surface gravity determinations for the primary 


co mponents and so we only consider the primaries from 
the Torres et al. (20101 sample. 

It IS important to note that the eclipsing binary sys¬ 
tems used for the surface gravity calibration are more dis¬ 
tant than the stars for which we can interferometrically 
determine angular diameters and effective temperatures 
for. Thus, for the surface gravity calibration it was nec¬ 
essary to compute the dereddened indices {b — y)o, toq, cq 
in order to obtain the highest accuracy possible for the 
intermediate-group stars, which rely on oq (an index us¬ 
ing dereddened colors) as a temperature indicator. No¬ 
tably, however, we found that the dereddened photome¬ 
try actually worsened \ogg determinations for the early 
and late groups. Dereddened colors were computed using 
the IDL routine UVBYBETA. 

The results of the logg calibration are presented in 
Table and Figure As described above, for the late 
group stars. (Teff <8500 K), logg is determined in the 
(b—y)—ci plane. The mean and median of the log resid¬ 
uals (in the sense of grid-fundamental) are - 0.001 dex and 
-0.0 38 d ex, respectively, and the RMS error 0.145 dex. As 
in § |4.H we found that the j3 — ci plane produced less ac¬ 
curate atmospheric parameters, relative to fundamental 
determinations, for late group stars. 

For the intermediate group stars (8500 K < Teff < 
11000 K), log 5 is determined in the oq — r* plane. The 
mean and median of the log 5 residuals are -0.060 dex 
and -0.069 dex, respectively with RMS error 0.091 dex. 
For the early group stars {T^s > 11000 K), log 5 is de¬ 
termined in the Ci — j3 plane. The mean and median of 
the log 5 residuals are -0.0215 dex and 0.024 dex, respec¬ 
tively, with RMS error 0.113 dex. The [u — b\~ fi plane 
was also investigated for early group stars, but was found 
to produce log g values of lower accuracy relative to the 
fundamental determinations. 

When considered collectively, the mean and median of 
the logg residuals for all stars are -0.017 dex and -0.034 
dex, and the RMS error 0.127 dex. The uncertainties 
in our surface gravities that arise from propagating the 
photometric errors through our atmospheric parameter 
determination routines are of the order ~ 0.02 dex, sig¬ 
nificantly lower than the uncertainties demonstrated by 
the comparison to fundamental values of logg. 

As stated above, the main concern with using double- 
lined eclipsing binaries as surface gravity calibrators for 
our photometric technique is contamination from the un¬ 
resolved secondary components. The logg residuals were 
examined as a function of both mass ratio and orbital 
period. While the amplitude of the scatter is marginally 
larger for low mass ratio or short period systems, in all 
cases our log 5 determinations are within 0.2 dex of the 
fundamental values « 85% of the time. 

To assess any potential systematic inaccuracies of the 
grids themselves, the surface gravity residuals were ex¬ 
amined as a function of Teff and the grid-determined 
log 5 . Figures show the logg residuals as a function 
of Teff and log 5 , respectively. No considerable system¬ 
atic effects as a function of either effective temperature 
or logg were found in the uvbyj5 determinations of log g. 

The effect of rotational velocity on our logg determi¬ 
nations was considered. As before, v sin i data for the 
surface gravity calib rators was collected from [Glebocki] 


& Gnacinski (20051). As seen in Figure 10 the major- 
ity of the log g calibrators are somewhat slowly rotating 
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(tisini < 150 km s“^). While the usini correction in¬ 
creases the accuracy of our logg determinations for the 
early-group stars in most cases, the correction appears 
to worsen our determinations for the intermediate group, 
which appear systematically high to begin with. 

The potential systematic effect of metallicity on our 
log 5 determinations is considered in Figure |11[ show¬ 
ing the surface gravity residuals as a function mpe/H]. 
Metallicity measurements were available for very few of 
these s t ars, and were primari l y take n from |Ammons et al.| 
(2006); Anderson & Francis (2012). Nevertheless, there 
does not appear to be a global systematic trend in the 
surface gravity residuals with metallicity. There is a 
larger scatter in log g determinations for the more metal- 
rich, late-type stars, however it is not clear that this effect 
is strictly due to metallicity. 

In summary, for the open cluster tests we assign logg 
uncertainties in three regimes: ± 0.145 dex for stars be¬ 
longing to the late group, ± 0.091 dex for the intermedi¬ 
ate group, and ± 0.113 dex for the early group. 

For our sample of nearby field stars we opt to assign 
a uniform systematic uncertainty of ± 0.116 dex for all 
stars. We do not attempt to correct for any systematic 
effects by applying offsets in logg, as we did with Tgg. 
As noted in discussion of the T^s calibration, we do ap¬ 
ply the V sin i correction to both intermediate and early 
group stars, as these corrections permit us to better re¬ 
produce open cluster ages (as presented in §|^. 

4.2.2. Comparison with Spectroscopic Measurements 

The Balmer lines are a sensitive surface gravity indi¬ 
cator for stars hotter than ^ 9000 K and can be 
used as a semi-fundamental surface gravity calibration 
for the early- and intermediate-group stars. The rea¬ 
son why surface gravities derived using this method are 
considered semi-fundamental and not fundamental is be¬ 
cause the method still relies on model atmospheres for 
fitting the observed line profiles. Nevertheless, surface 
gravities determined through this method are considered 
of high fidelity and so we performed an additional consis¬ 
tency check, comparing our uvbyl3 values of log g to those 
with well-determined spectroscopic logg measurements. 

N9 3 fit theoretica l profiles of hydrogen Balmer lines 
from Kurucz (|1979 ) to high resolution spectrograms of 
the H/t and JI 7 lines for a sample of 16 stars with 
uvbyP photometry. The sample of 16 stars was mostly 
drawn from the list of photometric /3 standards of |Craw-] 
ford] (|1966|) . We compared the log g values we deter¬ 
mined through interpolation in the uvbyjS color grids 
to the semi-fundamental spectroscopic values determined 
by N93. The results of this comparison are presented in 
Table [3 

Though N93 provide dereddened photometry for the 
spectroscopic sample, we found using the raw HM98 
photometry produced significantly better results (yield¬ 
ing an RMS error that was three times lower). For the 
early group stars, the atmospheric parameters were de¬ 
termined in both the cq — /3 plane and the [u — b] — jS 
plane. In both cases, /3 is the gravity indicator, but we 
found that the log 5 values calculated when using cq as 
a temperature indicator for hot stars better matched the 
semi-fundamental spectroscopic log 5 values. This result 
is consistent with the result from the effective temper¬ 
ature calibration that suggests Cq better predicted the 


effective temperatures of hot stars than [u — b]. As be¬ 
fore, logg for intermediate group stars is determined in 
the ao — r* plane. 

We tested uvbyP color grids of different metallicity, 
alpha-enhancement, and microturbulent velocity and de¬ 
termined that the non-alpha-enhanced, solar metallicity 
grids with microturbulent velocity i^turb = 0 km s“^ best 
reproduced the spectroscopic surface gravities for the 
sample of 16 early- and intermediate-group stars mea¬ 
sured by N93. 

The log (7 residuals, in the sense of (spectroscopic - 
grid), as a function of the grid-calculated effective tem¬ 
peratures are plotted in Figure There is no evidence 
for a significant systematic offset in the residuals as a 
function of either the uvbyjS-deteTmmed T^g or log 5 . For 
the early group, the mean and median surface gravity 
residuals are -0.007 dex and 0.004 dex, respectively, with 
RMS 0.041 dex. For the intermediate group, the mean 
and median surface gravity residuals are -0.053 dex and 
-0.047 dex, respectively, with RMS 0.081 dex. Consider¬ 
ing both early and intermediate group stars collectively, 
the mean and median surface gravity residuals are -0.027 
dex and -0.021 dex, and the RMS 0.062 dex. 

One issue that may cause statistically larger errors in 
the log g determinations compared to the T^g determina¬ 
tions is the linear interpolation in a low resolution log¬ 
arithmic space (the uvbyjS colors are calculated at steps 
of 0.5 dex in log 5 ). In order to mitigate this effect one 
requires either more finely gridded models or an inter¬ 
polation scheme that takes the logarithmic gridding into 
account. 

4.3. Summary of Atmospheric Parameter 
Uncertainties 

Precise and accurate stellar ages are the ultimate goal 
of this work. The accuracy of our ages is determined 
by both the accuracy with which we can determine at¬ 
mospheric parameters and any systematic uncertainties 
associated with the stellar evolutionary models and our 
assumptions in applying them. The precision, on the 
other hand, is determined almost entirely by the preci¬ 
sion with which we determine atmospheric parameters 
and, because there are some practical limits to how well 
we may ever determine T^g and logg, the location of the 
star in the H-R diagram (e.g. stars closer to the main 
sequence will always have more imprecise ages using this 
method). 

It is thus important to provide a detailed accounting of 
the uncertainties involved in our atmospheric parameter 
determinations, as the final uncertainties quoted in our 
ages will arise purely from the values of the , o’log g 
used in our calculations. Below we consider the con¬ 
tribution of the systematics already discussed, as well as 
the contributions from errors in interpolation, photome¬ 
try, metallicity, extinction, rotational velocity, multiplic¬ 
ity, and spectral peculiarity. 

Systematics: The dominant source of uncertainty in 
our atmospheric parame ter d eterm inati ons are the sys¬ 
tematics quantified in § |4.1| and § |4 2| All systematic 
effects inherent to the uvby/3 method, and the particu¬ 
lar model color grids chosen, which we will call tTsys, are 
embedded in the comparisons to the stars with funda¬ 
mentally or semi-fundamentally determined parameters, 
summarized as approximately ^ 3.1% in Tf.g and ^ 0.116 
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Fig. 8 .— Compari son of the uvhvQ der ived log 5 values with fundamental values for the primary components of the double lined eclipsing 
binaries compiled in|Torres et al.|i2010||. Red, teal, and black points represent late, intermediate, and early group stars, respectively. In 
each case the solid colored line represents the mean of the residuals, Alog^ (in the sense of fundamental-As can be seen, the mean 
offsets for the late and early groups is negligible. For the intermediate group, however, while only five stars were used for calibration, the 
uvhyp log 5 - values are about 0.13 dex lower than the fundamental values on average. 
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Fig. 9. — Surface gravity residuals, Alog^ (in the sense of 
fundamental-rtubt//?), as a function of un 6 y^-determined log(Teff) 
(left) and lo go fright). Sol i d poin ts represent eclipsing binary pri¬ 
maries from |Torres et al.| ||2010|l and open circles are stars with 
spectroscopic logg deterrninations in N93. Of the 39 eclipsing bi¬ 
naries, only six have residuals greater than 0.2 dex in magnitude. 
This implies that the uvby(3 grids determine logg to within 0.2 dex 
of fundamental values ~ 85% of the time. Surface gravity resid¬ 
uals are largest for the cooler stars. Photometric surface gravity 
measurements are in better agreement with spectroscopic determi¬ 
nations than the eclipsing binary sample. There is no indication 
for a global systematic offset in ut;i>y/3-determined logg values as 
a function of either Teff or log^. 
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Fig. 11. — Surface gravity residuals, Alog^ (in the sense of 
fundamental-un6^/5), as a function o f [Fe/Hl. The m e tallici ty val¬ 
ues have been taken primaril y from I Ammons et al . ( 200 61, with 
additional values coming from |Anderson 6c t^'rancisl ( 2U12||. While 
metallicities seem to exist for very tew ot the surface gravity cali¬ 
brators used here, there does not appear to be a systematic trend 
in the residuals with [Fe/H]. There is a larger amount of scatter for 
the more metal-rich late-type stars, however the scatter is confined 
to a relatively small range in [Fe/H] and it is not clear that this 
effect is due to metallicity effects. 
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Fig. 10. — Surface gravity residuals, Alogp- (in the sense of 
fundamental-wn6^/3), of eclipsing binary primaries as a function 
of V si n 2 . Arrows indicate the locations of points after application 
of the |Figueras Blasi| ( Il998| l i;sin 2 correction, where in this case 
late group stars received tne same correction as the intermediate 
group. 


dex in logg. We also found that for stars with avail¬ 
able [Fe/H] measurements, the accuracy with which we 
can determine atmospheric parameters using uvby(3 pho¬ 
tometry does not vary systematically with metallicity, 
though we further address metallicity issues both below 
and in an Appendix. 

Interpolation Precision: To estimate the errors in 
atmospheric parameters due to the numerical precision 
of the interpolation procedures employed here, we gen¬ 
erated 1000 random points in each of the three relevant 
uvbyP planes. For each point, we obtained ten indepen¬ 
dent Teff, logg determinations to test the repeatability of 
the interpolation routine. The scatter in independent de¬ 
terminations of the atmospheric parameters were found 
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TABLE 2 

Primary Components of Double-lined Eclipsing Binaries with fundamental determinations of log^. 


Star 

Sp. Type 

Tsfif 

(K) 

'^uvby 

(K) 

logSEB 

(dex) 

9uvby 

(dex) 

V sin i 
(km s“^) 

[Fe/H] 

(dex) 

(b-y) 

(mag) 

mi 

(mag) 

Cl 

(mag) 

d 

(mag) 

EM Car 

08V 

34000 ± 2000 

21987 

3.855 

± 

0.016 

3.878 

146.0 

— 

0.279 

-0.042 

0.083 

2.617 

V1034 Sco 

09V 

33200 ± 900 

28228 

3.923 

± 

0.008 

3.969 

159.0 

-1.0 

0.190 

-0.024 

-0.068 

2.587 

AH Cep 

BO.SVn 

29900 ± 1000 

24867 

4.017 

± 

0.009 

4.115 

154.0 

— 

0.290 

-0.064 

0.003 

2.611 

V578 Mon 

BIV 

30000 ± 740 

25122 

4.176 

± 

0.015 

4.200 

107.0 

— 

0.206 

-0.024 

-0.003 

2.613 

V453 Cyg 

B0.4IV 

27800 ± 400 

24496 

3.725 

± 

0.006 

3.742 

130.0 

— 

0.212 

-0.004 

-0.004 

2.590 

CW Cep 

B0.5V 

28300 ± 1000 

22707 

4.050 

± 

0.019 

3.716 

120.0 

— 

0.355 

-0.077 

0.050 

2.601 

V539 Ara 

B3V 

18100 ± 500 

17537 

3.924 

± 

0.016 

3.964 

85.6 

— 

-0.033 

0.089 

0.268 

2.665 

CV Vel 

B2.5V 

18100 ± 500 

17424 

3.999 

± 

0.008 

3.891 

42.8 

— 

-0.057 

0.083 

0.273 

2.659 

AG Per 

B3.4V 

18200 ± 800 

15905 

4.213 

± 

0.020 

4.311 

92.6 

-0.04 

0.048 

0.079 

0.346 

2.708 

U Oph 

B5V 

16440 ± 250 

15161 

4.076 

± 

0.004 

3.954 

350.0 

— 

0.081 

0.050 

0.404 

2.695 

V760 Sco 

B4V 

16900 ± 500 

15318 

4.176 

± 

0.019 

4.061 

— 

— 

0.169 

0.023 

0.392 

2.701 

GG Lup 

B7V 

14750 ± 450 

13735 

4.298 

± 

0.009 

4.271 

123.0 

— 

-0.049 

0.115 

0.514 

2.747 

C Phe 

B6V 

14400 ± 800 

13348 

4.121 

± 

0.004 

4.153 

111.0 

— 

-0.039 

0.118 

0.559 

2.747 

Hya 

B8V 

11750 ± 190 

11382 

3.710 

± 

0.007 

3.738 

131.0 

— 

-0.020 

0.110 

0.841 

2.769 

V906 Sco 

B9V 

10400 ± 500 

10592 

3.656 

± 

0.012 

3.719 

81.3 

— 

0.039 

0.101 

0.996 

2.805 

TZ Men 

AOV 

10400 ± 500 

10679 

4.224 

± 

0.009 

4.169 

14.4 

— 

0.000 

0.142 

0.918 

2.850 

V1031 Ori 

A6V 

7850 ± 500 

8184 

3.559 

± 

0.007 

3.793 

96.0 

— 

0.076 

0.174 

1.106 

2.848 

0 Aur 

Aim 

9350 ± 200 

9167 

3.930 

± 

0.005 

3.894 

33.2 

-0.11 

0.017 

0.173 

1.091 

2.889 

V364 Lac 

A4m: 

8250 ± 150 

7901 

3.766 

± 

0.005 

3.707 

— 

— 

0.107 

0.168 

1.061 

2.875 

V624 Her 

A3m 

8150 ± 150 

7902 

3.832 

± 

0.014 

3.794 

38.0 

— 

0.111 

0.230 

1.025 

2.870 

V1647 Sgr 

AIV 

9600 ± 300 

9142 

4.252 

± 

0.008 

4.087 

— 

— 

0.040 

0.174 

1.020 

2.899 

VV Pyx 

AIV 

9500 ± 200 

9560 

4.087 

± 

0.008 

4.004 

22.1 

— 

0.028 

0.161 

1.013 

2.881 

KW Hya 

A5m 

8000 ± 200 

8053 

4.078 

± 

0.006 

4.390 

16.6 

— 

0.122 

0.232 

0.832 

2.827 

WW Aur 

A5m 

7960 ± 420 

8401 

4.161 

± 

0.005 

4.286 

35.8 

— 

0.081 

0.231 

0.944 

2.862 

V392 Car 

A2V 

8850 ± 200 

10263 

4.296 

± 

0.011 

4.211 

163.0 

— 

0.097 

0.108 

1.019 

2.889 

RS Cha 

A8V 

8050 ± 200 

7833 

4.046 

± 

0.022 

4.150 

30.0 

— 

0.136 

0.186 

0.866 

2.791 

MY Cyg 

FOm 

7050 ± 200 

7054 

3.994 

± 

0.019 

3.882 

— 

— 

0.219 

0.226 

0.709 

2.756 

El Cep 

F3V 

6750 ± 100 

6928 

3.763 

± 

0.014 

3.904 

16.2 

0.27 

0.234 

0.199 

0.658 

2.712 

FS Mon 

F2V 

6715 ± 100 

6677 

4.026 

± 

0.005 

3.992 

40.0 

0.07 

0.266 

0.148 

0.594 

2.688 

PV Pup 

A8V 

6920 ± 300 

7327 

4.255 

± 

0.009 

4.386 

66.4 

— 

0.200 

0.169 

0.636 

2.722 

HD 71636 

F2V 

6950 ± 140 

6615 

4.226 

± 

0.014 

4.104 

13.5 

0.15 

0.278 

0.157 

0.496 

— 

RZ Cha 

F5V 

6450 ± 150 

6326 

3.905 

± 

0.006 

3.808 

— 

0.02 

0.312 

0.155 

0.482 

— 

BW Aqr 

F7V 

6350 ± 100 

6217 

3.979 

± 

0.018 

3.877 

— 

— 

0.328 

0.165 

0.432 

2.650 

V570 Per 

F3V 

6842 ± 50 

6371 

4.234 

± 

0.019 

3.998 

44.9 

0.06 

0.308 

0.165 

0.441 

— 

CD Tau 

F6V 

6200 ± 50 

6325 

4.087 

± 

0.007 

3.973 

18.9 

0.19 

0.314 

0.178 

0.436 

— 

V1143 Cyg 

F5V 

6450 ± 100 

6492 

4.322 

± 

0.015 

4.155 

19.8 

0.22 

0.294 

0.165 

0.451 

2.663 

VZ Hya 

F3V 

6645 ± 150 

6199 

4.305 

± 

0.003 

4.182 

— 

-0.22 

0.333 

0.145 

0.370 

2.629 

V505 Per 

F5V 

6510 ± 50 

6569 

4.323 

± 

0.016 

4.325 

31.4 

-0.03 

0.287 

0.142 

0.435 

2.654 

HS Hya 

F4V 

6500 ± 50 

6585 

4.326 

± 

0.005 

4.471 

23.3 

0.14 

0.287 

0.160 

0.397 

2.648 


Spectral type, temperature, and fundamental l og q information originate from [Torres et al. (|2010il. The uvbv0\oEQ val u es are from this 
work. Projected rotational velocities are from |Glebocki & Gnacinski| (|2005|l, [f 'e/tl] trom Siderson & Francis| (|2012||; |Ammons et al.| 
(|2006|l, and the uvbyp photometry are from HMyy. T he surtace gravities m the column log Quvby ^re derived together with 'I'uvby, and 
not for the Teff values given by [Torres et al.| | |2010[ |. 


to be < 10“^° K, dex, and thus numerical errors are 
assumed zero. 

Photometric Errors: Considering the most basic el¬ 
ement of our approach, there are uncertainties due to 
the propagation of photometric errors through our atmo¬ 
spheric parameter determination pipeline. As discussed 
in §[^ the photometric errors are generally small 0.005 
mag m a given index). Translating the model grid points 
in the rectangular regions defined by the magnitude of 
the mean photometric error in a given index, and then 
interpolating to find the associated atmospheric param¬ 
eters of the perturbed point, we take the maximum and 
minimum values for Teg and logg to calculate the error 
due to photometric measurement error. 

To simplify the propagation of photometric errors for 
individual stars, we performed simulations with ran¬ 
domly generated data to ascertain the mean uncertainty 
in Teff, logg that results from typical errors in each of 
the uvbyP indices. 

We begin with the HM98 photometry and associated 
measurement errors for our sample (3499 stars within 
100 pc, B0-F5, luminosity classes IV-V). Since the HM98 


compilation does not provide Oq or r*, as these quanti¬ 
ties are calculated from the four fundamental indices, we 
calculate the uncertainties in these parameters using the 
crude approximation that none of the uvbyP indices are 
correlated. Under this assumption, the uncertainties as¬ 
sociated with ao and r* are as follows: 

CTao = 1.362cr2_^ + -f (17) 

CT.* = ^0.07^al_y + 0.35^al+aj. (18) 

A model for the empirical probability distribution 
function (hereafter PDF) for the error in a given uvbyfi 
index is created through a normalized histogram with 25 
bins. From this empirical PDF, one can randomly draw 
values for the error in a given index. For each uvbyjS 
plane, 1,000 random points in the appropriate range of 
parameter space were generated with photometric errors 
drawn as described above. The eight (Teff, logg) values 
corresponding to the corners and midpoints of the “stan¬ 
dard error rectangle” centered on the original random 
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TABLE 3 

Stars with semi-fundamental determinations of logg through Balmer-line fitting. 


HR 

Sp. Type 

Tfr 

(K) 

'^uvby 

(K) 

9 spec 

(dex) 

Quvhy 

(dex) 

(.b-y) 

(mag) 

mi 

(mag) 

Cl 

(mag) 

[u- 6] 
(mag) 

/3 

(mag) 

63 

A2V 

8970 

9047 

3.73 

3.912 

0.026 

0.181 

1.050 

1.425 

2.881 

153 

B2IV 

20930 

20635 

3.78 

3.872 

-0.090 

0.087 

0.134 

0.264 

2.627 

1641 

B3V 

16890 

16528 

4.07 

4.044 

-0.085 

0.104 

0.319 

0.485 

2.683 

2421 

AOIV 

9180 

9226 

3.49 

3.537 

0.007 

0.149 

1.186 

1.487 

2.865 

4119 

B6V 

14570 

14116 

4.18 

4.176 

-0.062 

0.111 

0.481 

0.673 

2.730 

4554 

AOVe 

9360 

9398 

3.82 

3.863 

0.006 

0.155 

1.112 

1.425 

2.885 

5191 

B3V 

17320 

16797 

4.28 

4.292 

-0.080 

0.106 

0.297 

0.470 

2.694 

6588 

B3IV 

17480 

17025 

3.82 

3.864 

-0.065 

0.079 

0.292 

0.418 

2.661 

7001 

AOVa 

9540 

9508 

4.01 

3.977 

0.003 

0.157 

1.088 

1.403 

2.903 

7447 

B5III 

13520 

13265 

3.73 

3.712 

-0.016 

0.088 

0.575 

0.743 

2.707 

7906 

B9IV 

10950 

10838 

3.85 

3.861 

-0.019 

0.125 

0.889 

1.130 

2.796 

8585 

AIV 

9530 

9615 

4.11 

4.175 

0.002 

0.170 

1.032 

1.373 

2.908 

8634 

B8V 

11330 

11247 

3.69 

3.672 

-0.035 

0.113 

0.868 

1.077 

2.768 

8781 

B9V 

9810 

9868 

3.54 

3.593 

-0.011 

0.128 

1.129 

1.380 

2.838 

8965 

B8V 

11850 

11721 

3.47 

3.422 

-0.031 

0.100 

0.784 

0.969 

2.725 

8976 

B9IVn 

11310 

11263 

4.23 

4.260 

-0.035 

0.131 

0.831 

1.076 

2.833 


Spectral type, Teff, and spectroscopic log^ originate from N93. The uvby^ Teff and logp- values are from this work. Though N93 does 
not provide formal errors on the atmospheric parameters, those authors estimate uncertainties of ~ 0.03 dex in their spectroscopically 
determined log 5 -. The fractional errors in their photometrically derived Teff range from 2.5% for stars cooler than 5 ^ 11000 K to 4% for 
stars hotter than 20000 K. The photometry is from HM98. 


data point are then evaluated. The maximally discrepant 
(Teff, log 5 ) values are saved and the overall distributions 
of ATefi/Teff and A log g are then analyzed to assess the 
mean uncertainties in the atmospheric parameters de¬ 
rived in a given uvbyj5 plane due to the propagation of 
typical photometric errors. 

For the late group, points were generated in the range 
of {b — y) — Cl parameter space bounded by 6500 K 
< Teff < 9000 K and 3.0 < logg < 5.0. In this group, 
typical photometric uncertainties of (ub-y) = 0.003 mag 
and (ctci) = 0.005 mag lead to average uncertainties of 
0.6 % in Teff and 0.055 dex in logg. For the intermediate 
group, points were generated in the range of oq — r* pa¬ 
rameter space bounded by 8500 K < Teff < 11000 K and 
3.0 < logg < 5.0. In this group, typical photometric un¬ 
certainties of (cTao) = 0.005 mag and (cTr*) = 0.005 mag 
lead to average uncertainties of 0.8 % in Teff and 0.046 
dex in logg. For the early group, points were generated 
in the range of ci — /? parameter space bounded by 10000 
K < Teff < 30000 K and 3.0 < logg < 5.0. In this group, 
typical photometric uncertainties of {(7^) = 0.005 mag 
and (ap) — 0.004 mag lead to average uncertainties of 
1.1 % in Teff and 0.078 dex in logg. Across all three 
groups, the mean uncertainty due to photometric errors 
is « 0.9% in Teff and « 0.060 dex in logg. 

Metallicity Effects: For simplicity and homogene¬ 
ity, our method assumes solar composition throughout. 
However, our sample can more accurately be represented 
as a Gaussian centered at -0.109 dex with a « 0.201 
dex. Metallicity is a small, but non-negligible, effect 
and allowing [M/H] to change by ± 0.5 dex can lead 
to differences in the assumed Teff of ~ 1 - 2 % for late-, 
intermediate-, and some early-group stars, or differences 
of up to 6% for stars hotter than ~ 17000 K (of which 
there are few in our sample). In logg, shifts of ±0.5 dex 
in [M/H] can lead to differences of ^ 0.1 dex in the as¬ 
sumed log g for late- or early-group stars, or ~ 0.05 dex in 
the narrow region occupied by intermediate-group stars. 

Here, we estimate the uncertainty the metallicity ap¬ 
proximation introduces to the fundamental stellar pa¬ 


rameters derived in this work. We begin with the actual 
uvbyj3 data for our sa mple, and [Fe/H] measur ements 
from the XHIP catalog (Anderson & Francis 2012), which 
exist for approximately 68 % of our sample. Those au¬ 
thors collected photometric and spectroscopic metallicity 
determinations of Hipparcos stars from a large number 
of sources , calibrated the values to the high-resolution 
catalog of Wu et al. (2011 1 in an attempt to homogenize 
the various databases, ana published weighted means for 
each st ar. The calibration p r ocess is described in detail 
in §5 of Anderson & Francis (2012). 

For each ot the stars with available [Fe/H] in our field 
star sample, we derive Te®, log g in the appropriate uvbyfi 
plane for the eight cases of [M/H]=-2.5,-2.0,-1.5,-1.0,- 
0.5,0.0,0.2,0.5. Then, given the measured [Fe/H], and 
making the approximation that [M/H] = [Fe/H], we per¬ 
form a linear interpolation to find the most accurate val¬ 
ues of Teff,logg given the color grids available. We also 
store the atmospheric parameters a given star would be 
assigned assuming [M/H]=0.0. Figure [T^ shows the his¬ 
tograms of Teff/Teff jm/h]=o and logg-Tog g[M/H]=o- We 
take the standard deviations in these distributions to re¬ 
flect the typical error introduced by the solar metallicity 
approximation. For Teff, there is a 0.8% uncertainty in¬ 
troduced by the true dispersion of metallicities in our 
sample, and for logg, the uncertainty is 0.06 dex. These 
uncertainties in the atmospheric parameters are natu¬ 
rally propagated into uncertainties in the age and mass 
of a star through the likelihood calculations outlined in 

Reddening Effects: For the program stars studied 
here, interstellar reddening is assumed negligible^ _ Per¬ 

forming the reddening corrections (described in § | 2 . 2 [ ) on 
our presumably unreddened sample of stars within 100 
pc, we find for the ^ 80% of stars for which deredden¬ 
ing proved possible, that the distribution of Ay values in 
our sample is approximately Gaussian with a mean and 
standard deviation of g = 0.007, ct = 0.125 mag, respec¬ 
tively (see Figure [T9|). Of course, negative Ay values 
are unphysical, but applying the reddening corrections 
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Fig. 12.— Distributions of the true variations in Tgfr (left) and 
logg (right) caused by our assumption of solar metallicity. The 
“true” Tefi and logg values are determined for the ~ 68% of our 
field star sample with [Fe/H] measurements in XHIP and from 
linear interpolation between the se t of atmospheric parameters de- 
termined in eight ATLAS9 grids l|Castelli & Kurucz||2006l |2004|l 
that vary from -2.5 to 0.5 dex in [iVi/tlJ. ^ 



Fig. 13.— The effect of interstellar reddening on atmospheric 
parameters derived from uvbyfi photometry. The iso chrones and 
mass tracks plotted are those of|Bressan et al.| (|2012|l. The tail of 
each vector represents a given point m a specihc photometric plane 
{{b — y) —Cl for the late group stars in red, ag — r* for the intermedi¬ 
ate group stars in teal, and ci —/? for the early group stars in black) 
and its corresponding value in [Tefi, log g]. The tip of the vector 
points to the new value of [Tgjf, log g] after each point in photomet¬ 
ric space has been “dereddened” assuming arbitrary values of Ay ■ 
The shifts in uvbyfS space have been computed acco rding to the 
extinction measurem ents of|Schlegel et al.|||1998|l and|Crawford &| 
|Mandwewala| (|1976|l, assuming Ay ~ 4.ii3Yii/f6 — y). me magni¬ 
tudes ot Ay chosen for this figure represent the extremes of values 
expected for our sample of nearby stars and are meant to illustrate 
the directionality of the effects of reddening as propagated through 
the uvbyf} planes. Finally, note for the early group (black vectors), 
the Ay values are an order of magnitude larger and much higher 
than expected for our sample. Again, this is to illustrate the direc¬ 
tionality of the reddening effect, which is particularly small for the 
early group which rely on ci, the Balmer discontinuity index, for 
temperature, and /3, a color between two narrow-band filters with 
nearly the same central wavelength, for logg. 


to our uvbyj3 photometry and deriving the atmospheric 
parameters for each star in both the corrected and un¬ 
corrected cases gives us an estimate of the uncertainties 
in those parameters due to our assumption of negligible 
reddening out to 100 pc. The resulting distributions of 
Teff,o/2eff and log go “logg, where the naught subscripts 
indicate the dereddened values, are sharply peaked at 1 
and 0, respectively. The FWHM of these distributions 
indicate an uncertainty < 0.2% in Teff and ~ 0.004 dex 
in logg. For the general case of sources at larger dis¬ 
tances that may suffer more significant reddening, the 
systematic effects of under-correcting for extinction are 


illustrated in Figure [13 

Uncertainties iriProjected Rotational Veloci¬ 


ties: The Glebocki & Gnacinski| (2005) compilation con- 
tains mean usinz measurements, as well as individual 
measurements from multiple authors. Of the 3499 stars 
in our sub-sample of the HM98 catalog, 2547 stars have 
V sin i values based on 4893 individual v sin i measure¬ 
ments, 1849 of which have an accompanying measure¬ 
ment error. Of these measurements, 646 are for interme¬ 
diate or early groups, for which rotation corrections are 
performed in our method. The mean fractional error in 
vsini for this subset of measurements is ^ 13%. Caclu- 
lating the atmospheric parameters for these stars, then 
performing the FB98 v sin i corrections using ripot and 
Vrot ± ^v^at allows us to estimate the magnitude of the 
uncertainty in Teff, log g due to the uncertainties in v sin i 
measurements. The resulting RMS errors in Teff, log g are 
0.7% and 0.01 dex, respectively. When usini measure¬ 
ments are not available, an average value based on the 
spectral type can be assumed, resulting in a somewhat 
larger error. The systematic effects of under-correcting 
for rotation are illustrated in Figured 

Influence of Multiplicity: In a large study such as 
this one, a high fraction of stars are binaries or higher 
multiples. Slightly more than 30% of our sample stars 
are known as members of multiple systems. We choose 
not to treat these stars differently, given the unknown 
multiplicity status of much of the sample, and caution 
our readers to use due care regarding this issue. 

Influence of Spectral peculiarities: Finally, early- 
type stars possess several peculiar subclasses (e.g. Ap, 
Bp, Am, etc. stars) for which anomalous behavior has 
been reported in the uvbyj3 system with respect to their 
“normal-type” counterparts. Some of these peculiarities 
have been linked to rotation, which we do account for. 
We note that peculiar subclasses constitute ~ 4% of our 
sample and these stars could suffer unquantified errors 
in the determination of fundamental parameters when 
employing a broad methodology based on calibrations 
derived from mostly normal-type stars (see Tables 
for a complete accounting of the spectral types usM for 
calibrations). As these subclasses were included in the 
atmospheric parameter validation stage (§^, and sat¬ 
isfactory accuracies were still obtained, we mose not to 
adjust our approach for these stars and estimate the un¬ 
certainties introduced by their inclusion is negligible. 

Final Assessment: Our final atmospheric parame¬ 
ter uncertainti es ar e dom inat ed by the systematic effects 
quantified in § |4.1| and § |4.2[ with the additional effects 
outlined above contributing very little to the total un¬ 
certainty. The largest additional contributor comes from 
the photometric error. Adding in quadrature the sources 
f^sys, f^num, Tphot, sin i and cfAy results in final 

error estimates of 3.4% in Tes and 0.14 dex in logg. 

The use of uvbyj3 photometry to determine fundamen¬ 
tal stellar parameters is estimated in previous literature 
to lead t o uncertainties of j ust 2.5% in Teg and 0.1 dex 
in logg (Asiain et al. 1997), with our assessment of the 
errors somewhat higher. 

The uncertainties that we derive in our Stromgren 
method work can be compared with those given by 
other methods. The Geneva photometry system 
(t7, Rl, R2, Ul, G filters), like the Stromgren system, has 
been used to derive TpfF.logg, and [M/H] values based 
on atmospheric grids (Kobi fc North|[l990 Kunzli et al. 
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Fig. 14.— Comparison of PARSEC isochrones (solid lines), Ek- 
strom isochrones in the rotating case (dashed lines), and Ekstrom 
isochrones in the non-rotating case (dotted lines). The solid black 
lines are evolutionary tracks for stars of intermediate-mass, from 
the PARSEC models. All evolutionary tracks plotted are for solar 
metallicity. 


19971, with Kunzli et al. (1997) finding 150-250 K (few 
percent) errors in log'Jeff and 0. 1-0.15 dex errors in log 
g, comparable to our values. From stellar model at¬ 
mosphere fitting to high dispersion spectra, errors of 1 - 
5% in Teff and 0.05-0.15 dex (typically 0.1 dex) in logq 
are quoted for early type stars (e.g. |Nieva & Simon- 
Di'az| 2011|) , though systematic effects in log g on the 
order ot an additional 0.1 dex may be present. Wu et 
al. ( 2011 ) tabulate the dispersions in atmospheric pa¬ 
rameters among many different studies, finding author- 
to-author values that differ for OBA stars by 300-5000 
K in Teff (3-12%) and 0.2-0 .6 dex in logg (cm/s^), and 
for FGK stars 40-100 K in Tgff and 0.1-0.3 dex in logg 
(cm/s^). 

5. AGE ESTIMATION FROM ISOGHRONES 
5.1. Selection of Evolutionary Models 

Once Teff and log 5 have been established, ages are 
determined through a Bayesian grid search of the fun¬ 
damental parameter space encompassed by the evolu¬ 
tionary models. In this section we discuss the selection 
of evolution models, the Bayesian approach, numerical 
methods, and resulting age/mass uncertainties. 

Two sets of isochrones are considered in this work. The 
model families are compared i n Figurefldl The PA RSEG 
solar-metallicity isochrones of|Bressan et al. (2012), here¬ 
after B12, take into account in a sell-consistent manner 
the pre-main-sequence phase of evolution. The PARSEC 
models are the most recent iteration of the Padova evo¬ 
lutionary models, with significant revisions to the major 
input physics such as the equation of state, opacities, 
nuclear reaction rates and networks, and the inclusion of 
microscopic diffusion. The models are also based o n the 
new reference so lar composition, Z — 0.01524 from Caf- 
fau et al. ( 2011 ), but can be generated for a wide range 
of metallicities. The B12 models cover the mass range 
0.1 - 12 M©. 

PARSEC isochrones are attractive because early-type 
dwarfs have relatively rapid evolution with the pre-main- 
sequence evolution constituting a significant fraction of 
their lifetimes, i.e. tpms/tms is larger compared to 


stars of later types. For stars with effective temper¬ 
atures in the range 6500 K - 25000 K (approximately 
spectral types B0-F5), the B12 models predict pre-main 
sequence lifetimes ranging from ^ 0.2-40 Myr, main- 
sequence lifetimes from ^14 Myr - 2.2 Gyr, and the ra¬ 
tio rpMs/Tvis ~ 1-6 — 2.4%. A star of given initial mass 
thus can be followed consistently through the pre-MS, 
MS, and post-MS evolutionary stages. As a consequence, 
most points in Teff —log g space will have both pre-ZAMS 
and post-ZAMS ages as possible solutions. Figure [^il¬ 
lustrates the evolution of atmospheric and corresponding 
photometric properties according to the PARSEC mod¬ 
els. 


The solar-metallicity isochrones of Ekstrom et al. 
([2012 ), hereafter E12, also use updated opacities and nu¬ 


clear reaction rates, and are the first to take into account 
the effects of rotation on global stellar properties at inter¬ 
mediate masses. They are available for both non-rotating 
stars and stars that commence their lives on the ZAMS 
with a rotational velocity of 40% their cri tical rotational 


velocit y (urot,i/ccrit = 0.4); however, the Ekstrom et al. 
( 2012 ) models do not take the pre-main sequence phase 


into account. The E12 models currently exist only for so¬ 
lar metallicity (Z=0.014 is used), but cover a wider range 
of masses ( 0.8 — I 2 OM 0 ). 

The E12 models are attractive because they explictly 
account for rotation, though at a fixed percentage of 
breakup velocity. All output of stellar evolutionary mod¬ 
els (e.g. lifetimes, evolution scenarios, nucleosynthesis) 
are affected by axial stellar rotation which for massive 
stars enhances the MS lifetime by about 30% a nd may in¬ 
crease isochro nal age estimates by about 25% (jMeynet &| 
Maeder|2000). In terms of atmospheres, for A-type stars. 


stellar rotation increases the strength of the Balmer dis¬ 
continuity re lative to a non-rotating star w ith the same 
color index (Maeder & Peytremann 1970|. In the El 2 
models, the convective overshoot parameter was selected 
to reproduce the observed main sequence width at inter¬ 
mediate masses, which is important for our aim of dis¬ 
tinguishing the ages of many field stars clustered on the 
main sequence with relati vely large uncertainties in their 
surface gravities. Figure [T^ shows, however, that there 
is close agreement betweentne B12 and the rotating E12 
models. Thus, there is not a significant difference be¬ 
tween the two models in regards to the predicted width 
of the MS band. 

It should be noted that the uvbyfd grids of |GasteiII| 


& Kurucz (2006 2004) were generated assuming a solar 
metallicity value ot Z=0.017. As discussed elsewhere, 
metallicity effects are not the dominant uncertainty in 
our methods and we are thus not concerned about the 
very small metallicty differences between the two model 
isochrone sets nor the third metallicity assumption in the 
model atmospheres. 

In matching data to evolutionary model grids, a gen¬ 
eral issue is that nearly any given point in an H-R di¬ 
agram (or equivalently in Teg-logg space), can be re¬ 
produced by multiple combinations of stellar age and 
mass. Bayesian inference can be used to determine the 
relative likelihoods of these combinations, incorporating 
prior knowledge about the distributions of the stellar pa¬ 
rameters being estimated. 
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5.2. Bayesian Age Estimation 

A simplistic method for determining the theoretical 
age and mass for a star on the Hertzsprung-Russell (H- 
R) diagram is interpolation between isochrones or evolu¬ 
tionary model s. Some problems w i th this approach, a s 
pointed out by|Takeda et al.|(|2007|);|Pont & Eyer|(|2004|), 
is that interpolation betweeii isochrones neither accounts 
for the nonlinear mapping of time onto the H-R diagram 
nor the non-uniform distribution of stellar masses ob¬ 
served in the galaxy. As a consequence, straightforward 
interpolation between isochrones results in an age distri¬ 
bution for held stars that is biased towards older ages 
compared to the distribution predicted by stellar evolu¬ 
tionary theory. 

Bayesian inference of stellar age and mass aims to 
eliminate such a bias by accounting for observationally 
and/or theoretically motivated distribution functions for 
the physical parameters of interest. As an example, for a 
given point with error bars on the H-R diagram, a lower 
stellar mass should be considered more likely due to the 
initial mass function. Likewise, due to the longer main- 
sequence timescales for lower mass stars, a star that is 
observed to have evolved off the main sequence should 
have a probability distribution in mass that is skewed 
towards higher masses, i.e. because higher mass stars 
spend a more signihcant fraction of their entire lifetime 
in the post-MS stage. 


5.2.1. Bayes Formalism 

Bayesian estimation of the physical parameters can 
proceed from comparison of the data with a selection 
of models. Bayes’ Theorem states: 


P(model|data) ex P(data|model) x P(model) (19) 

The probability of a model given a set of data is pro¬ 
portional to the product of the probability of the data 
given the model and the probability of the model itself. 
In the language of Bayesian statistics, this is expressed 
as: 


posterior ex likelihood x prior. (20) 

Our model is the set of stellar parameters, age (r) and 
mass (M*), and our data are the measured effective tem¬ 
perature, Teff, and surface gravity, logg, for a given star. 
At any given combination of age and mass, the predicted 
Teff and log g are provided by stellar evolutionary models. 
The statistic for an individual model can be computed 
as follows: 


x’(r,M.) = 5:>k^ (21) 

_ [{Tes)o - (reff)i;]^ [(logg)o - (logg)i;]^ 

( 7 ^ (J? ’ 

Aff log g 

( 22 ) 

where the subscripts O and E refer to the observed and 
expected (or model) quantities, respectively, and a is the 
measurement error in the relevant quantity. 

Assuming Gaussian statistics, the relative likelihood of 
a specific combination of (Teff,logg) is: 


P(data|model) = P(reff,obs, loggobslr, M*) 


(X exp 


1 

2 


X^(t, M*) 


(23) 

(24) 


Finally, the joint posterior probability distribution for 
a model with age r and mass M*, is given by: 


P(model|data) = P(r, M*|Teff,obs, log gobs) (25) 


oc exp 


1 

2 


x"(t,M*) 


P(t)P(M*), 

(26) 


where P(t) and P(M*) are the prior probability dis¬ 
tributions in age and mass, respectively. The prior prob¬ 
abilities of age and mass are assumed to be independent 
such that P(t, M*) = P(r)P(M*). 


5.2.2. Age and Mass Prior Probability Distribution 
Funetions 

Standard practice in the Bayesian estimation of stellar 
ages is to assume an ag e prior that is uniform in linear 
age, e . g. |Pont fc Ey er (|2004|); [jprgensen &: Li ndegren 
(|2005|);|i’akeda et al. 'PD(j7p;|Melsen et al.|(|2U13p. 'i'here 
are two main justihcations tor choosing a uniform age 
prior: 1) it is the least restrictive choice of prior and 2) at 
this stage the assumption is consistent with observations 
that suggest a fairly constant star f ormation rate in sola r 
neighborhood over the past 2 Gyr (Cignoni et al.|2006). 


Since the evolutionary models are logarithmically grid- 
ded in age, the relative probability of age bin i is given 
by the bin width in linear age divided by the total range 
in linear age: 


^’(log(T,) < log(T) < l0g(Ti+i)) = 


'^n '^O 


(27) 


where r„ and tq are the largest and smallest allowed 
ages, respectively. This weighting scheme gives a uniform 
probability dist ribution in lin e ar age . 

As noted by |Takeda et al.| (|2007| , it is important to 
understand that assuming a flat prior in linear age cor¬ 
responds to a highly non-uniform prior in the measured 
quantities of log Teg and logg. This is due to the non¬ 
linear mapping between these measurable quantities and 
the physical quantities of mass and age in evolutionary 
models. Indeed, the ability of the Bayesian approach to 
implicitly account for this effect is considered one of its 
main strengths. 

As is standard in the Bayesian estimation of stellar 
masses, an initial mass function (IMF) is assumed for 
the prior probability distribution of all possible stellar 
masses. Several authors point out that Bayesian esti¬ 
mates of physical parameters are relatively insensitive 
to the mass prior (i.e. the precise form of the IMF as¬ 
sumed), especially in the case of parameter determina¬ 
tion over a small or moderate range in mass space. For 
this work cons idering BAF stars, the power law IMF of 
Salpeter (1955) is assumed for the mass prior, so that the 
relative probability of mass bin i is given by the following 
expression: 










































18 


-2.35 


P(M, <M < M,+i) oc Ml 
5.2.3. Numerical Methods 


(28) 


As Takeda et al. (2007) point out, in Bayesian age es- 
timation interpolation should be performed only along 
isochrones and not between them. To avoid biasing our 
derived physical parameters from interpolating between 
isochrones, we generated a dense grid of PARSEC mod¬ 
els. The evolutionary models were acquired with a spac¬ 
ing of 0.0125 dex in log(age/yr) and 0.0001 Mq in mass. 
All probabilities were then computed on a 321 x 321 grid 
ranging from log(age/yr )=6 to 10 and from 1-10 Mq. 

5.2.4. Age and Mass Uncertainties 

Confidence intervals in age and mass are determined 
from the one-dimensional marginalized posterior prob¬ 
ability distributions for each parameter. Since the 
marginalized probability distributions can often be as- 
symetric, the region chosen for determining confidence 
intervals is that of the Highest Posterior Density (HPD). 
This method selects the smallest range in a parame¬ 
ter that encompasses N% of the probability. The HPD 
method is discussed in more detail in the appendix. 

Notably, uncertainties in the ages depend on where 
in the log 5 and logTgff parameter space the star is lo¬ 
cated, and whether a pre-main sequence or a post-zero- 
age-main sequence age is more appropriate. In the pre- 
main sequence phase, both atmospheric parameters are 
important in age determination. For post-ZAMS stars, 
however, the relative importance of the two parameters 
changes. When stars are just bouncing back from the 
ZAMS and are starting to evolve through the MS phase, 
log 5 must be known precisely (within the range of ~4.3 
to 4.45) in order to derive a good age estimate. The age 
at which this bounce occurs will be a function of mass 
(earlier for more massive stars). Otherwise, once late B, 
A, and early F stars are comfortably settled on the MS, 
their e volu tion is at roughly constant temperature (see 
Figure 14) and so the gravity precision becomes far less 


important, with temperature precision now critical. 

6 . THE METHODOLOGY TESTED ON OPEN 
CLUSTERS 

An important test of our methods is to assess the ages 
derived from our combination of uvbyj3 photometry, at¬ 
mospheric parameter placement, and comparison to evo¬ 
lutionary models relative to the accepted ages for mem¬ 
bers of well-studied open clusters. We investigate four 
such clusters with rigorous age assessment in previous lit¬ 
erature: IC 2602, a Persei, the Pleiades, and the Hyades. 

The youngest (< 20 — 30 Myr) open clusters may be 
age-dated kinematically, by tracing the space motions 
of individual members back to the time whe n the stars 
were in closest proximity to one another (Soderblom 
Mot . After < 1 galactic rotation period, however, m- 
dividual member motions are randomized to the extent 
of limiting the utility of the kinematic method. Be¬ 
yond ~ 20 — 30 Myr, the most precise open cluster ages 
come from the lithium depletion boundary (LDB) tech¬ 
nique. This method uses the lithium abundances, which 
diminish predictably with time, of the lowest mass clus¬ 
ter members to converge on precise (^10%) ages. LDB 



Fig. 15.— Histograms of the visual extinction, in mag¬ 

nitudes for individual members of the four open clusters consid¬ 
ered here. The extinction values are calculated using the relation 
Ay = A.237E (b — y), with the (b — y) color excesses computed as 
described in § |2.2| 


ages are available for IC 2602: 
et al . 2010), a Per: r = 90 ± 10 Myr (I Stauffer et 


T = 461^ Myr (Dobbie 


1999 1, and the Pleiades: r = 125 ±8 Myr (jbtautter et al. 


19981. The LDB technique does not work past ~ 250 
Myr, so the Hyades is dated based on isochrone fitting in 
the H-R diagram using stars with high precision distance 
measurements, with cu rrently accepted age 625 ±50 Myr 
( Perryman et al.||1998 ). 


6.1. Process 

Membership probabilities, uvbyfd photometry, and pro¬ 
jected rotational velocities are obtained for members 
of these open clusters via the WEBDA open cluster 
databas^ For the Pleiades, membershi p information 
was au gmented and cross-referenced with [Stauffer et al.| 
(|2007|) . Both individual uvbyj3 measurements and calcu- 
lations of the mean and scatter from the literature mea¬ 
surements are available from WEBDA in each of the pho¬ 
tometric indices. As the methodology requires accurate 
classification of the stars according to regions of the H-R 
diagram, we inspected the spectral types and /3 indices 
and considered only spectral types B0-F5 and luminosity 
classes HI-V for our open cluster tests. 

In contrast to the field stars studied in the next sec¬ 
tion, the open clusters studied here are distant enough for 
interstellar reddening to significantly affect the derived 
stellar parameters. The photometry is thus dereddened 
as described in § |2.2[ Figure [TS] shows the histograms of 
the visual extinction Ay for each cluster, with the impact 
of extinction on the atmospher ic p arameter determina¬ 
tion illustrated above in Figure [T^ 

In many cases, individual cluster stars have multiple 
measurements of wsinz in the WEBDA database and 
we select the measurement from whichever reference is 
the most inclusive of early-type members. In very few 
cases does a cluster member have no rotational veloc¬ 
ity measurement present in the database; for these stars 
we assume the mean v sin i accordin g to the Tf.fr — v sin i 
relation presented in Appendix B of Gray| (2005). 

Atmospheric parameters are determined for each clus¬ 
ter member, as described in § [^ Adopting our knowl¬ 
edge from the comparison to fundamental and semi- 

® http://www.miivie.ac.at/webda/ 
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Fig. 16.— PARSEC isochrones and mass tracks (|Bressan et al. 120121) in logTgff-logg space and the isochrones of |Ekstrdm et aL](|2012|| 
(including rotation, plotted as dashed lines) with our uvbyj3 pnotometric determinations of the atmospheric parameters, i-or early and 
intermediate group stars, the black filled circles represent the usini corrected atmospheric parameters (using the FB98 formulae), while the 
open circles represent the uncorrected parameters. Note that the late-group stars do not receive a vsini correction but are still plotted as 
filled circles. In both cases the point sizes are oc usini. The typical uncertainties in our logTeff and logg determinations are represented by 
the error bars at the bottom of the figure. These uncertainties correspond to 1.6% or ^ 0.007 dex in logTeff and 0.091 dex (intermediate), 
and 0.145 dex (late) in log 5 -, corresponding to the RMS errors as determined in the ef fective temperature and surface gravity calibrations. 
Top left: IC 2602 members; the currently accepted age of IC 2602 is r = 46^*]® Myr l Dobbie et al.||2010 1 . Top right: Members of the a 
Persei cluster, which has a currently accepted age of r = 90 di 10 Myr dS t auffer et ai.||1999D. Bottom left: Pleiades members where the 
currently accepted age of the Pleiades is r = 125 di 8 Myr|Stauffer et al.| (|iyyy||. Ut the ~ 20 Pleiads that sit below the zero age main 
sequence, 5 are known pulsators of the 5 Sen or 7 Dor variety. Additionally, there is an excess of slow rotators sitting below the ZAMS. 
Possible reasons for this observed behavior include systematics of the atmospheric models (several authors have noted problems with the 
treatment of convection in ATLAS9 models at this mass range), failure of the evolutionary models to predict the true width of the main 
sequence (though this effect is unlikely to be as large as the scatter seen here), and overaggressi ve dereddening procedu res. Bottom right: 
Hyades cluster members where the currently accepted age of the H yades is r = 625 di 5 0 Myr iPerryman et ah ||1998|l. Note the far left 
outlier, HD 27962, is a known blue straggler (|Abt|1985| |Eggen|1995|l and was ex cluded by I Perryman et al.| (TujiiJa l m tneir isochrone-fitting 
analysis. The outlier far below the ZAMS, HiJ 27268, is a spectroscopic binary | |Debernarcli et al.|2UUu| r 
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fundamental atmospheric parameters (§ 4.1 & § 4.2), a 
uniform 1.6% shift towards cooler Teg was applied to all 
temperatures derived from the model color grids to ac¬ 
count for systematic effects in those grids. The FB98 
V sin i corrections were then applied to the atmospheric 
parameters. The wsinz corrections prove to be a cru¬ 
cial step in achieving accurate ages for the open clusters 
(particularly for the Pleiades). 


6.2. Results 

The results of applying our procedures to open clus¬ 
ter samples appear in Figure While the exact 

cause(s) of the remaining scatter observed in the empiri¬ 
cal isochrones for each cluster is not known, possible con¬ 
tributors may be systematic or astrophysical in nature, 
or due to incorrect membership information. Multiplic¬ 
ity, variability, and spectral peculiarities were among the 
causes investigated for this scatter, but the exclusion of 
objects on the basis of these criteria did not improve age 
estimation for any individual cluster. The number of 
stars falling below the theoretical ZAMS, particulary for 
stars with log Teg < 3.9, is possibly systematic and may 
be due to an incomplete treatment of convection by the 
ATLAS9 models. T his so urce of uncertainty is discussed 
in further detail in § |8.2| 

For each cluster, we publish the individual stars con- 
sid ered , al ong with relevant parameters, in Tables i i 
[H & El In each table, the spectral types and v sin i 
measurements are from WEBDA, while the dereddened 
uvbyP photometry and atmospheric parameters are from 
this work. 


6.2.1. Ages from Bayesian Inferenee 

Once atmospheric parameters have been determined, 
age determination proceeds as outlined in § [^ For each 
individual cluster member, the likelihood) and pos¬ 
terior probability distribution are calculated for each 
point on a grid ranging from log(age/yr)=6.5 to 10, with 
masses restricted to 1 < M/Mq < 10. The resolu¬ 
tion of the grid is 0.0175 dex in log(age/yr) and 0.045 
Mq in mass. The 1-D marginalized posterior PDFs for 
each individual cluster member are normalized and then 
summed to obtain an overall posterior PDF in age for the 
cluster as a whole. This composite posterior PDF is also 
normalized prior to the determination of statistical mea¬ 
sures (mean, median, confidence intervals). Additionally, 
the posterior PDFs in log(age) for each member are mul¬ 
tiplied to obtain the total probability in each log(age) bin 
that all members have a single age. While the summed 
PDF depicts better the behavior of individual stars or 
groups of stars, the multiplied PDF is best for assigning 
a single age to the cluster and evaluating any potential 
systematics of the isochrones themselves. 

As shown in Figure the summed age PDFs for 
each cluster generally follow the same behavior: (1) the 
peaks are largely determined by the early group (B-type) 
stars which have well-defined ages due to their unambigu¬ 
ous locations in the T^ff — \ogg diagram; (2) examining 
the age posteriors for individual stars, the intermediate 
group stars tend to overpredict the cluster age relative 
to the early group stars, and the same is true for the 
late group stars with respect to the intermediate group 
stars, resulting in a large tail at older ages for each of 


the summed PDFs due to the relatively numerous and 
broad PDFs of the later group stars. For IC 2602 and 
the Pleiades, the multiplied PDFs have median ages and 
uncertainties that are in close agreement with the liter¬ 
ature ages. Notably, the results of the open cluster tests 
favor an age for the Hyades that is older (~ 800 Myr) 
than the accepted value, though not q uite as old as the 
recent estimate of 950±100 Myr from Brandt fc Huang] 
(2015). The Bayesian age analysis also favors an age for 
a Per that is younger (~ 70 Myr) than the accepted value 
based on lithium depletion, but older than the canonical 
50 Myr fro m the Upper Main Sequence Turnoff [Mermil-| 
hod (1981). In an appendix, we perform the same anal- 
ysis tor the open clusters on p(t) rather than p(logT), 
yielding similar results. 

The results of the open cluster test are presented in 
Table H It is noted that all statistical measures of the 
marginalized age PDFs quoted hereafter are from PDFs 
normalized in log(age), as opposed to converting to lin¬ 
ear age and then normalizing. This choice was made 
due to the facts that 1) the isochrones are provided in 
uniform logarithmic age bins, and 2) the marginalized 
PDFs of individual stars are more symmetric (and thus 
better characterized by traditional statistical measures) 
in log(age) than in linear age. Notably, the median age is 
equivalent regardless of whether one chooses to analyze 
prob(logT) or prob(T). This issue is discussed further in 
an appendix. In general, there is very close agreement 
in the Bayesian method ages between B12 and rotating 
E12 models. For IC 2602 and the Pleiades, our analy¬ 
sis yields median cluster ages (as determined from the 
multiplied PDFs) that are within 1-cr of accepted values, 
regardless of the evolutionary models considered. The 
Bayesian analysis performed with the PARSEC models 
favor an age for a Persei that is ^ 20% younger than the 
currently accepted value, or ^ 20% older for the Hyades. 

6.2.2. Ages from Isochrone Fitting 

As a final test of the two sets of evolutionary mod¬ 
els, we used -minimization to find the best-fitting 
isochrone for each cluster. By fitting all members of a 
cluster simultaneously, we are able to assign a single age 
to all stars, test the accuracy of the isochrones for stel¬ 
lar ensembles, and test the ability of our uvby/3 method 
to reproduce the shapes of coeval stellar populations in 
Teff—log (? space. For this exercise, we did not interpolate 
between isochrones, choosing instead to use the default 
spacing for each set of models (0.1 dex and 0.0125 dex in 
log(age/yr) for the E12 and B12 models, respectively). 
For the best results, we consider only the sections of the 
isochrones with logg between 3.5 and 5. 0 de x. The re¬ 
sults of this exercise are shown in Figure 18 The best¬ 
fitting El 2 isochrone (including rotation) is consistent 
with accepted ages to within 1% for the Pleiades and 
Hyades, ^ 15% for IC 2602, and ^ 44% for a-Per. For 
the B12 models, the best-fit isochrones are consistent 
with accepted ages to ^ 8% for the Pleiades, ~ 20% for 
the Hyades and IC 2602, and ~ 47% for a-Per. The 
B12 models produce systematically younger ages than 
the E12 models, by a fractional amount that increases 
with absolute age. 

As detailed above, the open cluster tests revealed that 
our method is able to distinguish between ensembles of 
differing ages, from tens to hundreds of Myr, at least 
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Fig. 17.— Left panels: ID marginalized, normalized posterior PDFs in age, calculated from|Bressan et aL|||2012|| evolutionary models, for 
individual open cluster members. Black, teal, and red histograms represent early, intermediate, and late group stars, respectively. Middle 
panels: Sums of the individual PDFs depicted on the left. This figure shows the total probability associated with the 200 age bins between 
log(age/yr)=6.5 to 10. The grey shaded regions indicate the currently accepted ages of IC 2602 (46^® Myr), a Per (OOdilO Myr), the 
Pleiades (125di8 Myr), and the Hyades (625di50 Myr). Right panels: Products of the individual PDFs depicted in the left panels. The 
grey shaded regions again depict the accepted literature age ranges of each cluster. 


TABLE 4 

Open Cluster Ages 


Cluster 

Lit. Age 
(Myr) 

Models 


Summed PDF 
Median 
(Myr) 

Summed PDF 
68 % C.I. 
(Myr) 

Multiplied PDF 
Median 
(Myr) 

Multiplied PDF 
68 % C.I. 

(Myr) 

Xmin 

(Myr) 

IC 2602 

46t® 

Ekstrom et al. 

(2012 

1 80 

32-344 

42 

41-46 

39 



IBressan et al.l 

120121 

79 

27-284 

46 

44-50 

37 

a Persei 

Ekstrom et al. 

(2012 

1 234 

83-1618 

71 

68-74 

50 


125^3 

liriressan et al.l 

I2U12I 

226 

74-1500 

70 

69-74 

48 

Pleiades 

Ekstrom et al. 

(2012 

1 277 

81-899 

128 

126-130 

126 



lusressan et al.l 

|2U12| 

271 

85-948 

123 

121-126 

115 

Hyades 

Ekstrom et al. 

(2012 

1 872 

518-1940 

827 

812-837 

631 



[Bressan et al.| 

2012 

844 

487-1804 

764 

747-780 

501 


Literature ages (column 2) come from the sources referenced in §[m For each set of evolutionary models, the median and 68% confidence 
interval are computed for both the summed PDF (columns 4,5) ana multiplied PDF (columns 6,7). The final column indicates the best-fit 
isochrone found through x^-mi^iimization of all cluster members in log (Teff) ~ Mgp space. Note, the Hyades analysis includes the blue 
straggler HD 27962 and the spectroscopic binary HD 27268. Excluding these outliers results in a median and 68% confidence interval of 
871 Myr [517-1839 Myr] of the summed PDE or 832 Myr [812-871 Myr] of the multiplied PDE, using the B12 models. 
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in a statistical sense. For individual stars, large uncer¬ 
tainties may remain, particularly for the later types, ow¬ 
ing almost entirely to the difficulty in determining both 
precise and accurate surface gravities. The open cluster 
tests also demonstrate the importance of a u sin i correc¬ 
tion for early (BO-AO) and intermediate (A0-A3) group 
stars in determining accurate stellar parameters. While 
the usinj correction was not applied to the late group 
(A3-F5 in this case) stars, it is likely that stars in this 
group experience non-negligible gravity darkening. The 
typically unknown inclination angle, i, also contributes 
significant uncertainties in derived stellar parameters and 
hence ages. 


7. THE METHODOLOGY APPLIED TO NEARBY 
FIELD STARS 

As an application of our developed, calibrated, val¬ 
idated, and tested methodology, we consider the com¬ 
plete HM98 photometric catalog of 63,313 stars. We are 
interested only in nearby stars that are potential tar¬ 
gets for high contrast imaging campaigns, and for which 
interstellar extinction is negligible. We thus perform a 
distance cut at 100 pc, using dista nces from the XHIP 
catalog (Anderson & Francis 2012). We perform an ad¬ 
ditional cut in spectral type (using information from 
XHIP), considering only B0-F5 stars belonging to lumi¬ 
nosity classes IV,V, because this is the range for which 
our method has been shown to work with high fidelity 
and additionally these are the primary stars of interest 
to near-term high-contrast imaging surveys. In total, we 
are left with 3499 stars. Figure shows the distribu¬ 
tion of our field star sample in spectral type, distance. 
Ay, [Fe/H], and wsinz. The distributions of photomet¬ 
ric errors in given uvby/3 indices are shown in Figure 
^ 0 l and the mean errors in each index are summarized 
as follows: {(Jb-y) , {(Jmi) , , {^ 0 ) , {<^ao) = 

0.003,0.004,0.005,0.004,0.005,0.005 mag. 

Projected rotational velocities f or the sample of nearby 
field s tars are sourced from the IGlebocki fc Gnacinskil 
(2005) compilation, which contains v sin i measurements 
for 2874 of the stars, or ~ 82% of the sample. For an 
additional 8 stars v sin i measurements are collected from 


Zorec & Royer (2 012), and for another 5 stars ?;sinj 
values come from Schroder et al. (2009). For the re¬ 


maining stars without vsmi measurements, a projected 
rotational velocity is assumed according to th e mean 
V sin i — Teff relation from Appendix B of Gray (2005). 
Atmospheric parameters are corrected for rotational ve- 
locity effects as outlined in § |3.3| 

Atmospheric parameter determination was not possi¬ 
ble for six stars, due to discrepant positions in the rele¬ 
vant uvby/3 planes: HIP 8016 (a B9V Algol-type eclips¬ 
ing binary), HIP 12887 (a poorly studied F3V star), HIP 
36850 (a well-studied AlV-|-A2Vm double star system), 
HIP 85792 (a well-studied Be star, spectral type B2Vne), 
HIP 97962 (a moderately studied B9V star), and HIP 
109745 (an AOHl star, classified in XHIP as an AlIV 
star). Gonsequently, ages and masses were not computed 
for these stars. 

An H-R diagram of the entire sample is shown in Fig- 
ure m with the evolutionary models of |Bressan et al.| 
(2012|) overlaid. Equipped with atmospheric parameters 
lor the remaining 3493 stars, and assuming uniform un¬ 
certainties of 3.4% and 0.14 dex in and logfjf, respec¬ 


tively, ages and masses were computed via the process 
outlined in § [^ Posterior probab ilities were calculated 
on a uniform 321x321 grid of thefBressan et al. (2012) 
models, gridded from 1 Myr-10 Gyr in steps ot (J.0125 
dex in log (age), and from 1 -lOMr:^ in steps of O.O 28 M 0 . 
As the Bressan et al. (2012) models exist for high resolu¬ 
tion timesteps, no interpolation between isochrones was 
required. 

From the 2D joint posterior PDF, we obtain the 
marginalized ID PDFs in age and mass, from which we 
compute the mean (expected value), median, mode (most 
probable value), as well as 68 % and 95% confidence in¬ 
tervals. Interpolated ages and masses are also included, 
and these values may be preferred, particularly for ob¬ 
jects with an interpolated age < 10 ® yr and a logg plac¬ 
ing it near the ZAMS (see § 8.2 for more detail). The 
table of ages and masses for all 3943 stars, including our 
newly derived atmospheric parameters, are available as 
an electronic table and a portion (sorted in ascending 
age) is presented here in Table In rare instances (for 
^ 5% of the sample), true 68 % and 95% confidence in¬ 
tervals were not obtained due to numerical precision, the 
star’s location near the edge of the computational grid, 
or some combination of the two effects. In these cases 
the actual confidence interval quoted is noted as a flag 
in the electronic table. 

As with the open clusters, we can sum the individual, 
normalized PDFs in age to produce composite PDFs for 
various subsets of our sample. Figure depicts the 
composite age PDF for our entire sample, as well as age 
PDFs for the subsets of B0-B9, A0-A4, A5-A9, and FO- 
F5 stars. From these PDFs we can ascertain the sta¬ 
tistical properties of these subsets of solar neighborhood 
stars, which are presented in Table 


7.1. Empirical Mass-Age Relation 

From our newly derived set of ages and masses of solar- 
neighborhood B0-F5 stars, we can determine an empiri¬ 
cal mass-age relation. Using the mean ages and masses 
for all stars in our sample, we performed a linear least 
squares fit using the NumPy polyfit routine, yielding the 
following relation, valid for stars 1.04 < M/Mq < 9.6: 


log(age/yr) = 9.532 - 2.929log • (29) 

The RMS error between the data and this relation is a 
fairly constant 0.225 dex as a function of stellar mass. 

7.2. Empirical Spectral-Type-Age/Mass Relations 

We can also derive empirical spectral-type-age and 
spectral-type-mass relations for the solar neighborhood, 
using the mean masses derived from our ID marginal¬ 
ized posterior PDFs in age, and spectral type informa¬ 
tion from XHIP. These relations are plotted in Figure 
|23[ and summarized in Table [7j 

8 . DISGUSSION 

The precision of the age-dating method described here 
relies on the use of Strbmgren ubvyj5 photometry to finely 
distinguish stellar atmosphere parameters and compare 
them to isochrones from stellar evolution models. For 
ages < 10 Myr and > 100 Myr, in particular, there is 
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Fig. 18.— Best fitting isochrones found through y^- minimization for four o pen clusters, with atmospheric parameters deter mined through 
uvhvd photometry. Left panels are the fully rotating |Ekstrdm et al.l ||2012|l evolutionary models while right panels are the [Bressan et al.| 
(|2012|| models . For the Pleiades, the best fitting isochrone age (i2t) Adyrj from the E12 models is within the currently accepted range ot 
l25±8 Myr. The B12 models give a best-fit age of 115 Myr, representing a fractional error of ~ 8% (or 1.25cr) relative to the accepted 
age. In the case of the Hyades (lower panels), the low and far left outliers are a spectroscopic binary and a blue straggler, respectively. 
Excluding these stars yields no change in the best-fitting isochrone for the E12 models and only moderately increases the best-fitting B12 
model to 530 Myr. 
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[Fe/H] 




spectral type 


Fig. 19. — Characterization of our sam ple of 3499 nearby field star s. Upper panels: histograms of the spectral types (left) and distances 
(right) of stars in our sample, taken from | Anderson S z Fra n(^ l|2012||. Middle panels: histogr ams of the V-band ex t inctio n in magnitudes 
(left), as derived by the IDL program described m ^ IMT and tne [i'e/H] values in dex fro m [Anderson &: Francis] l|201^ . Lower panels: 
histogram of the projected rotational velocities in our sample (left), with data taken from[Sebocki & Linacmski| ||2DU5|l, and usini as a 
function of spectral type (right) with grey x’s indicating individual stars and black squares representing the mean usmz in each spectral 
type bin. The error ba rs represent the standard deviation in v sini values for each bin. The red triangles indicate the empirical Teff-usini 
relation of jGrayj | |2005| | using the spectral-type-Teff relation o f jHabets &: Heintzej p'98T] | . 
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TABLE 5 

Ages, Masses, and Atmospheric Parameters of Nearby B0-F5 Field Stars 


HIP 

Tefi 

logs 

Mean 

Median 

Mode 

68 % 

95% 

Interp. 

Mean 

Median 

Mode 

68 % 

95% 

Interp. 




Age 

Age 

Age 

Age 

Age 

Age 

Mass 

Mass 

Mass 

Mass 

Mass 

Mass 


(K) 

(dex) 

(Myr) 

(Myr) 

(Myr) 

(Myr) 

(Myr) 

(Myr) 

(M©) 

(M©) 

(M©) 

(M©) 

(M©) 

(M©) 


65474 

24718 

4.00 

6 

7 

9 

5-12 

2-14 

13 

9.59 

9.61 

9.62 

9.4-9.9 

9.2-10.0 

10.26 

61585 

21790 

4.32 

6 

7 

11 

4-18 

1-21 

1 

7.53 

7.52 

7.48 

7.3-7.7 

7.1-8.0 

7.34 

61199 

16792 

4.18 

18 

22 

33 

13-53 

3-60 

36 

4.84 

4.83 

4.78 

4.6-5.0 

4.5-5.2 

4.95 

60718 

16605 

4.35 

19 

23 

35 

13-55 

3-61 

1 

4.75 

4.74 

4.70 

4.6-4.9 

4.4-5.1 

4.58 

60000 

15567 

4.12 

21 

26 

40 

14-65 

4-77 

60 

4.27 

4.26 

4.22 

4.1-4.4 

4.0-4.6 

4.48 

100751 

17711 

3.94 

23 

29 

43 

20-50 

5-52 

48 

5.41 

5.42 

5.35 

5.1-5.6 

5.0-5.9 

5.91 

23767 

16924 

4.10 

23 

30 

44 

18-56 

4-61 

46 

4.96 

4.95 

4.92 

4.7-5.2 

4.5-5.4 

5.14 

92855 

19192 

4.26 

24 

29 

34 

23-38 

8-40 

8 

6.39 

6.37 

6.25 

6 .0-6.6 

5.8-7.1 

5.95 

79992 

14947 

3.99 

26 

31 

48 

18-78 

4-89 

88 

4.01 

4.00 

3.97 

3.8-4.1 

3.7-4.3 

4.45 


The fractional uncertainty in our Tag determination s is 3.4% and the uncertainty in our logg determinations is 0.14 dex. All ages and 
masses are computed from the |Bressan et | |201^ models. Statistical measures are quoted for marginalized PDFs in log(age) rather 
than age, e.g. column 4 is 10^*°®*-rather than (t). Confidence intervals are computed via the HPD method. The full table containing 
ages, masses, and atmospheric parameters for all 3493 stars is available electronically. Table is published in its entirety in the electronic 
edition of ApJ, A portion is shown here for guidance regarding its form and content. 


TABLE 6 

Statistics of Composite Age PDFs 


Sp. Types 

Mean Age 
(Myr) 

Median Age 
(Myr) 

Mode Age 
(Myr) 

68 % C.I. 
(Myr) 

95% C.I. 
(Myr) 

B0-B9 

93 

122 

147 

56-316 

8-410 

A0-A4 

296 

365 

392 

200-794 

39-1090 

A5-A9 

572 

750 

854 

434-1372 

82-1884 

F0-F5 

1554 

1884 

2024 

1000-4217 

307-6879 




0.01 0.02 0.03 
CTb-y (mag) 


0.01 0.02 0.03 
ffao (mag) 




0.01 0.02 0.03 
o-r- (mag) 



(mag) 


(mag) 


ap (mag) 


Fig. 20.— Histograms of the uncertainties (in mag) for different 
uvbyp indices for the sample of ~ 3500 field stars discussed in 
§ The solid lines in each plot indicate the position of the mean 
uncertainty in that parameter. Uncertainties in ao and r* are 
calculated according to Eqns. (13) & (14). 


rapid evolution of log Teff and log g for intermediate-mass 
stars (see Figure [T]). This enables greater accuracy in age 
determination through isochrone placement for stars in 
this mass and age range. Fundamentally, our results rely 
on the accuracy of both the stellar evolution models and 
the stellar atmosphere models that we have adopted. Ac¬ 
curacy is further set by the precision of the photometry, 
the derived atmospheric parameters, the calibration of 
the isochrones, and the ability to determine whether an 
individual star is contracting onto the main sequence or 
expanding off of it. By using isochrones that include both 
pre- MS and post-MS evo lution in a self-consistent man¬ 
ner (Bressan et al. 2012), we can determine pre-ZAMS 


TABLE 7 

Empirical Spectral-Type Relations for Main Sequence 
B0-F5 Stars 


Sp. Type 

(r> 

(Myr) 

CTr 

(Myr) 

(M) 

(Me) 

(Me) 

No. of Stars 

BO 

19 

— 

4.75 

— 

1 

B1 

6 

— 

9.59 

— 

1 

B2 

15 

13 

6.96 

0.81 

2 

B3 

41 

16 

5.22 

0.50 

3 

B4 

26 

12 

4.94 

0.59 

4 

B5 

44 

16 

3.94 

0.49 

5 

B6 

84 

51 

3.69 

0.23 

4 

B7 

140 

209 

3.23 

0.60 

13 

B8 

99 

43 

3.28 

0.38 

18 

B9 

154 

86 

2.88 

0.88 

67 

AO 

285 

437 

2.47 

0.40 

120 

Al 

313 

217 

2.23 

0.29 

132 

A2 

373 

320 

2.11 

0.30 

144 

A3 

462 

412 

2.07 

0.96 

100 

A4 

540 

333 

1.84 

0.19 

37 

A5 

514 

350 

1.86 

0.81 

81 

A6 

628 

265 

1.85 

0.49 

46 

A7 

574 

262 

1.78 

0.30 

79 

A8 

642 

272 

1.64 

0.11 

62 

A9 

800 

339 

1.62 

0.21 

102 

FO 

994 

544 

1.52 

0.19 

324 

FI 

948 

352 

1.51 

0.13 

68 

F2 

1280 

526 

1.42 

0.19 

441 

F3 

1687 

633 

1.34 

0.23 

605 

F4 

1856 

600 

1.30 

0.12 

129 

F5 

2326 

697 

1.27 

0.18 

905 


in addition to post-ZAMS ages for every data point in 

Teff,logg). 

Above, we have described our methodology in detail, 
including corrections for reddening and rotation, and we 








































26 



5.0 

4.. 


• Sun 
^ p Pic 
^ (f And 
^ fDel 

★ Vega 

■ HR 8799 
^ Fomalhaut 


4.1 4.0 

log(Teff/K) 


Fig. 21.— H-R diagram for our sample of B0-F5 field stars within 100 pc. Thirteen stars with logg < 2.9 are excluded in this figure. 
Several stars of interest are plotted in gold. As before, red, teal, and black scatter points correspond to late, intermediate, and early group 
stars, respectively. Values for the Sun are also plotted for reference. Of note, ~ 770 of the stars plotted are subgiants according to their 
XHIP luminosity classes, while only ~ 250 stars have logg < 3.8, suggesting some spectral types are in error. 



Fig. 22. — Normalized composite age PDFs for our sample of field 
B0-F5 stars within 100 pc. The normalized composite PDFs are 
created by summing the normalized, ID marginalized age PDFs 
of individual stars in a given spectral type grouping. The black 
curve represents the composite pdf for all spectral types, while the 
colored curves represent the composite PDFs for the spectral type 
groups B0-B9, A0-A4, A5-A9, F0-F5 (see legend). Circles repre¬ 
sent the expectation values of the composite PDFs, while squares 
represent the medians. The solid and dashed lines represent the 
68% and 95% confidence intervals, respectively, of the composite 
PDFs. The statistical measures for these composite PDFs are also 
presented in Table 

have presented quality control tests that demonstrate the 
precision and accuracy of our ages. In the section we de¬ 
scribe several aspects of our analysis of specific interest, 
including the context of previous estimates of stellar ages 


for early type stars (§ 8.11, how to t reat stars with loca¬ 
tions apparently below the ZA MS 8.21, and discussion 
of notable individual objects (§ |8.3| . VVe will in the fu¬ 
ture apply our methods to new spectrophotometry. 


8.1. Methods Previously Employed in Age 
Determination for Early Type Stars 

In this section we place our work on nearby open clus¬ 
ter stars and approximately 3500 nearby held stars in 
the context of previous age estimation methods for BAF 
stars. 



isochrone placement in the H-R diagram, and literature 
ages (mostly coming from earlier application of uvbyP 
photometric methods). 

Among later type F dwarfs, previous age estimates 
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Fig. 23.— Empirical spectral-type-age relation (left) and spectral-type-mass relation (right) for solar neighborhood B0-F5 stars. Grey x’s 
represent individual stars, while the black scatter points represent the mean value in a given spectral type bin and the error bars represent 
the scatter in a that bin. 


come primarily from t he Geneva-Copenhagen Survey 
(Casagrande et al.|2011), but their reliability is caveated 
by the substantially diner ent values published in various 


iterations of the catalog ([Nordstrom et a h 200 4 [Holm- 
berg et al. 2007 2009 Casagrande et al. 2011) and the 


inherent dithculty ot .sochrone dating these later type 
dwarfs. 


More recently, Nielsen et al. (2013) applied a Bayesian 
inference approach to the age determination of 70 B- & 
A-type field stars via My vs B — V color-magnitude dia¬ 
gram isochrone placement, assuming a constant star for- 
mati on rate in the solar n eighborhood and a Salpeter 
IMF. De Rosa et al. (2014) estimated the ages of 316 A- 
type kars through placemient in a Mk vsV — K color- 
magnitude diagram. Both of these broad-band p hoto¬ 
metri c stud ies used the theoretical isochrones of ISiessI 
eTallpOOO). 

Considering the above sources of ages, the standard 
deviation among them suggests scatter among authors 
of only 15% for some stars up to 145%, with a typical 
value of 40%. The full range (as opposed to the disper¬ 
sion) of published ages is 3-300%, with a peak around 
the 80% level. The value of the age estimates presented 
here resides in the large sample of early type stars and 
the uniform methodology applied to them. 

8.2. Stars Below the Main Sequence 

In Figure [2T] it may be noted that many stars, particu¬ 
larly those with log Teff < 3.9, are located well below the 
model isochrones. Using rotation-corrected atmospheric 
parameters, ~ 540 stars or ~ 15% of the sample, fall 
below the theoretical ZAMS. 

Prior studies also faced a similarl y large fraction of 
stars falling below the main sequence. |Song et al.| (|2001[) 
arbitrarily assigned an age of 50 Myr to any star lying 
below the 10 0 Myr isochrone used in that work. |Tetzlaff| 
et al. (2011) arbitrarily shifted stars towards the ZAMS 
and treated them as ZAMS stars. 

Several possibilities might be invoked to explain the 
large population of stars below the log 5 — log Teff 
isochrones: ( 1 ) failure of evolutionary models to predict 
the true width of the MS, (2) spread of metallicities, with 
the metal-poor MS residing beneath the solar-metallicity 
MS, (3) overaggressive correction for rotational velocity 


effects, or (4) systematics involved in surface gravity or 
luminosity determinations. Of these explanations, we 
consider (4) the most likely, w i th (3) also contributing 
somewhat. Valenti & Fischer (2005) found a 0.4 dex 
spread in log 5 among their niain sequence FGK stars 
along with a 0.1 dex shift downward relative to the ex¬ 
pected zero metallicity main sequence. 

The Bayesian age estimates for stars below the MS are 
likely to be unrealistically old, so we compared the ages 
for these stars with interpolated ages. Using the field 
star atmospheric parameters and Bressan et al. (2012) 
models, we performed a 2D linear interpolation with 
the SciPy routine griddata. Stars below the main se¬ 
quence could be easily identified by selecting objects with 
log (age/yr)B,^y^,, - log (age/yr)i„jg,.p > 1.0. Notably, for 
these stars below the MS, the linear interpolation pro¬ 
duces more realistic ages than the Bayesian method. A 
comparison of the Bayesian and interpolated ages for all 
stars is presented in Figure Of note, there is closer 
agreement between the Bayesian and interpolation meth¬ 
ods in regards to estimating masses. 

Figure [^further serves to illustrate the difference be¬ 
tween the Bayesian ages and the interpolated ages, which 
scatters over an order magnitude from a 1:1 relationship. 
A number of stars that fall below the MS and have inde¬ 
pendently constrained ages are examined in detail in n 
§ 8.3 These stars have interpolated ages that are more 
in line with prior studies, and in light of this, we publish 
the interpolated ages in addition to the Bayesian ages in 
the final electronic table. 


8.3. Stars of Special Interest 

In this section we discuss stars of particular inter¬ 
est given that they have either spatially resolved debris 
disks, detected possibly pla netary mass compani ons, or 
both. As a final test of the Bressan et al. (2012) evolu¬ 
tionary models and our Bayesian age and mass estima¬ 
tion method, we performed our analysis on these stars, 
including the Sun. 


8.3.1. Sun 

The atmospheric parameters of our Sun are known 
with a precision that is orders of magnitude higher than 
what is obtainable for nearby field stars. One would thus 
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Fig. 24.— Comparison of ages for BAF field stars derived through 
2D linear interpolation and Bayesian inference. Grey points rep¬ 
resent those stars with Alogage/yr > 1 (in the sense of Bayesian 
minus interpolated), which coincide with the same stars that reside 
below the MS. 

expect the assumed priors to have a negligible influence 
on the Bayesian age and mass estimates. 

The effective temperature of the Sun is calculated to 
be Teff = 5771.8 ± 0.7 K from the total solar irradi- 
ance ( Kopp fc L ean 2011), the solar radius (Haberre- 
fne lAU 2009 definition of the AU, 


iter et al. 




constant. The solar surface gravity is calculated to be 
log 5 = 4.43812 ± 0.00013 dex from the lAU 2009 value 
of GMq and the solar radius (Haberreiter et al. 2008). 
Using these values, our Bayesian analysis yields a me¬ 
dian age of 5.209±0.015 Gyr. The Bayesian estimation 
also yields a mass estimate of 0.9691±0.0003 Mq. Per¬ 
forming a 2D linear interpolation yields a slightly older 
age of 5.216 Gyr and slightly lower mass of 0.9690 Mq. 
As expected, the precise solar values lead to an elliptical 
joint posterior PDF in age and mass, and symmetric ID 
marginalized PDFs. The difference between the Bayesian 
age estimate and interpolated age is negligible in this 
regime of extremely s mall uncertainti e s. Th is test also 
demonstrates that the Bressan et al. (2012) evolution¬ 


ary models may introduce a systematic overestimation 
of ages and underestimation of masses towards cooler 
temperatures, though because the Sun is substantially 
different from our sample stars we do not extrapolate 
this conclusion to our sample. 

8.3.2. HR 8199 

HR 8799 is located near th e ZAMS and is met al-poor 
with [Fe/H]=0.47± 0.10 dex (|Gray & Kaye|1999|). How¬ 
ever, because HR 8799 is a A Boo peculiar-type star, its 
photospheric metallicity may not reflect the global stel¬ 
lar metal abundance. The age of HR 8799 is believed to 
be 30li!^ Myr based on its proposed membership to the 
Columba association (Zuckerman et al.||2011 ). 

Figure [21] shows that HR 8799 lays well below the 
theoretical ZAMS. This location is well-documented from 
other spectroscopic and photometric analyses of the star, 
and is likely due to a combination of its genuine youth 
and subsolar metallicity. Gonsistent with the discussion 


in §8.2 and as illustrated in Figure [24] our Bayesian age 
analysis leads to an unrealistically old age for the star, 
with a median age of 956 Myr and a 68 % confidence 
interval of 708-1407 Myr. The Bayesian approach also 
seems to overestimate the mass, with a median mass of 
1 . 59 Af 0 and 68 % confidence interval of 1.49 — 1 . 68 M 0 . 
Notably, however, 2D linear interpolation leads to more 
reasonable age estimates: 26 Myr assuming our newly 
derived atmospheric parameters (Teff=7540 K, log 5 = 
4.43) , or 25 Myr using T eff=7430 K and logg = 4.35 
from [Gray &: Kaye (1999). 


8.3.3. /3 P%c 


Zuckerman (2001) assigned an age of 12 Myr to f3 Pic 
based on its proposed membership to the moving group 
of the same name. Isochronal age estimates f or the star 
have ranged from the Z AMS age to 300 Myr (Barrado y 
Navascues et al. 1999). Nielsen et al. (2013) performed 
a Bayesian analysis concluding a median age of 109 Myr 
with a 68 % confidence interval of 82-134 Myr. Although 
barely below the ZAMS, /3 Pic in our own Bayesian anal¬ 
ysis has a much older median age of 524 Myr with a 68 % 
confidence interval of 349-917 Myr. Prior authors also 
have noted that /3 Pic falls below the ZAMS on a color- 
magnitude diagram. As was the case for HR 8799, we 
conclude that our erroneous age for j3 Pic is due to the 
dominance of the prior assumption/s in exactly such a 
scenario. 

However, the interpolated age using our atmospheric 
parame ters of Tsff=8300 K, log (/=4.389, is 20 Myr. Us¬ 
ing the Gray et al.j (2006) values of Tefr=8052 K (within 
Icr of our determination), logg=4.15 (> 1.5cr away from 
our surface gravity) we obtain an interpolated age of 308 
Myr. 

8.3.4. K And 

K Andromedae is another proposed mem ber of the 
Columba association ( Zuckerman et a l.|201 1 ). Using the 
nominal 30 Myr age, [Carson et al.| (|2013|) suggested a 
companion discovered via direct ima g ing is of planetary 


mass (12 — 13Mjup). Hinkley et al. ([2013[) refuted this 
claim, concluding an age ot ZZU± 100 Myr from multiple 
isochronal analyses in §3.2 of that work. This older age 
estimate leads to a model-dependent companion mass of 
hO^iJgMjup. Our Bayesian analysis allows us to nearly 
rule out a 30 Myr age with a 95% confidence interval of 
29-237 Myr. The mean, median, mode, and 68 % confi¬ 
dence interval of the ID marginalized posterior PDF in 
age for k And are 118, 150, 191, and 106-224 Myr, re¬ 
spectively. Notably, k And has a projected rotational 
veloc ity of usinf ^ 160 km s“^ (Glebocki & Gnacinskif 
2005), and we find its rotation corrected atmospheric pa- 
rameters = 11903 ± 405 K, log g = 4.35 ± 0.14 dex) 
produce an interpolated age of 16 Myr. Using uncor¬ 
rected atmospheric parameters (Teff = 11263 ± 383 K, 
log 5 = 4.26 ± 0.14 dex) leads to an interpolated age of 
25 Myr. 

8.3.5. C Delphini 


De Rosa et al. (2014) recently published the discovery 
of a wide companion to C Delphini (HIP 101589). Those 
authors estimated the age of the system as 525±125 
Myr, from the star’s positions on a color-magnitude 
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log(age/yr) prob(mass) ^9® (Qy) prob(mass) 


Fig. 25.— 2D joint posterior PDFs in age and mass for an early-type star with typical atmospheric parameter uncertainties {left) and 
the Sun {right), for which Tgfj and logg are known to high precision. The dark, medium, and light blue shaded regions indicate the 68%, 
95%, and 99% confidence contours. Above, ID marginalized and normalized posterior PDF in age, with the shaded regions representing 
the same corresponding confidence intervals. Right, the same as above for mass. 


and a temperature-luminosity diagram, leading to a high-precision atmospheric parameters, 
model-dependent companion mass of 50±15 Mj^p. Our 
method yields a mean age of 552 Myr, with 68% and 
95% confidence intervals of 531-772 Myr, and 237-866 
Myr, respectively. Our rev ised age is in agreem ent 
with the previous estimate of De Rosa et al. (2014), al¬ 
though favoring the interpretation of an older system 
and thus more massive companion. The interpolated 
ages for C Del are somewhat older: 612 Myr for the 
rotation-corrected atmospheric parameters Tefj=8305 K, 
log 5=3.689, or 649 Myr for the uncorrected parameters 
Teff=8639 K, log 5=3.766. Note, in this case moderate 
rotation (usinf = 99.2 km s“^) leads to a discrepancy of 
only « 6 % in the derived ages. 


8.3.6. 49 Ceti 

49 Ceti does not have a known companion at present, 
but does po ssess a resolved molecular gas disk (Hughes 
et al. 2008). The star is a proposed member of the 40 
Myr Argus association, which would require cometary 
collisions to explain the gaseous disk that shoul d have 
dissipated by ^ 1 0 Myr due to radiation pressure (Zuck- 
erman et al. 120121) ■ With a mean rotational velocity of ^ 
190 km s“^ (Glebocki & Gnacinski||2005 ), and evidence 
that the star is highly inclined to our line of sight, ro¬ 
tational effects on photometric H-R diagram placement 
are prominent. Our uvbyj3 atmospheric parameters for 
49 Ceti are = 10007 ± 340 K, log 5 = 4.37 ± 0.14 
dex, after rotational effects were accounted for. These 
parameters place the star essentially on the ZAMS, with 
an interpolated age of 9 Myr, and calling into question 
the cometary genesis of its gaseous disk. However, the 
uncorrected atmospheric parameters (Teg = 9182 ± 309 
K, log 5 = 4.22 ± 0.14 dex) are more consistent with 
the Al spectral type and produces an interpolated age 
of 57 Myr, which seems to support the cometary colli¬ 
sion hypothesis. This case illustrates the importance of 


9. CONCLUSIONS 

In the absence of finely calibrated empirical age indica¬ 
tors, such as th e rotation-activity-age relation for solar- 
type stars (e.g. Mamajek & Hillenbrand 2008), ages for 
early spectral type stars typically have come from open 
cluster and moving group membership, or through as¬ 
sociation with a late-type companion that can be age- 
dated through one of the applicable empirical methods. 
Because of their rapid evolution, early type stars are 
amenable to age dating via isochrones. In this paper 
we have investigated the use of Stromgren photometric 
techniques for estimating stellar atmospheric parameters, 
which are then compared to isochrones from modern stel¬ 
lar evolution models. 

Bayesian inference is a particularly useful tool in the 
estimation of parameters such as age and mass from evo¬ 
lutionary models for large samples that span considerable 
ranges in temperature, luminosity, mass, and age. The 
Bayesian approach produces unbiased ages relative to 
a straightforward interpolation among isochrones which 
leads to age estimates that are biased towards older 
ages. However, as noted earlier, stars located beyond 
the range of the theory (below the theoretical ZAMS in 
our case) are assigned unreasonably old ages with the 
Bayesian method. This presumably is due to the clus¬ 
tering of isochrones and the dominance of the prior in in¬ 
ference scenarios in which the prior probability is chang¬ 
ing quickly relative to the magnitude of the uncertainty 
in the atmospheric parameters. Linear interpolation for 
stars apparently below the MS may produce more rea¬ 
sonable age estimates. 

The most important parameter for determining precise 
stellar ages near the ZAMS is the luminosity or surface 
gravity indicator. Effective temperatures, or observa¬ 
tional proxies for temperature, are currently estimated 
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with suitable precision. However, logg, luminosity, or 
absolute magnitude (requiring a precise distance as well), 
are not currently estimated with the precision needed to 
meaningfully constrain the ages of field stars near the 
ZAMS. This effect is particularly pronounced for lower 
temperature stars where, for a given shift in logg, the in¬ 
ferred age can change by many orders of magnitude. Our 
open cluster tests indicated that the age uncertainties 
due to the choice of evolutionary models are not signifi¬ 
cant compared to those introduced by the uncertainties 
in the surface gravities. 

We have derived new atmospheric parameters (taking 
stellar rotation into account) and model-dependent ages 
and masses for 3493 BAF stars within 100 pc of the 
Sun. Our method of atmospheric parameter determi¬ 
nation was calibrated and validated to stars with funda¬ 
mentally determined atmospheric parameters. We fur¬ 
ther tested and validated our method of age estimation 
using open clusters with well-known ages. In determining 
the uncertainties in all of our newly derived parameters 
we conservatively account for the effects of systematics, 
metallicity, numerical precision, reddening, photometric 
errors, and uncertainties in vsini as well as unknown 
rotational velocities. 

Field star ages must be considered with caution. At 
minimum, our homogeneously derived set of stellar ages 
provides a relative youth ordering. For those stars below 
the MS we encourage the use of interpolated ages rather 
than Bayesian ages, unless more precise atmospheric pa¬ 
rameters become available. Using the new set of ages, we 
presented an empirical mass-age relation for solar neigh¬ 
borhood B0-F5 stars. We also presented empirical rela¬ 
tions between spectral type and age/mass and we dis¬ 
cussed ages in detail for several famous low mass com¬ 
panion and/or debris disk objects. An anticipated use 
of our catalog is in the prioritization of targets for direct 
imaging of brown dwarf and planetary mass companions. 
David & Hillenbrand (2015b, in preparation) will explore 
how ages derived using this methodology can be applied 
to investigations such as debris disk evolution. 

The authors wish to thank John Stauffer for his help¬ 
ful input on sources of uvby/3 data for open clusters and 
Timothy Brandt for helpful discussions during the proof 
stage of this work regarding the open cluster analysis, re¬ 
sulting in the appendix material concerning logarithmic 
versus linear approaches and a modified version of Fig¬ 
ure 17. This material is based upon work supported in 
2014 and 2015 by the National Science Foundation Grad¬ 
uate Research Fellowship under Grant No. DGE1144469. 
This research has made use of the WEBDA database, op¬ 
erated at the Institute for Astronomy of the University of 
Vienna, as well as the SIMBAD database and VizieR cat¬ 
alog access tool, operated at GDS, Strasbourg, France. 
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APPENDIX 
Metallicity Effects 

We do not account explicitly for metallicity in this study, having assumed solar values in both our atmospheric models 
and our evolutionary grids. Our analysis in the Tgg and log 5 calibrations found that for stars with fundamentally 
determined atmospheric parameters and available [Fe/H] measurements, the accuracy with which we can determine 
atmospheric parameters using uvbyP photometry does not vary systematically with metallicity. _ 

The effects of different metallicity assumptions on the Stromgren index atmospheric grids i s illus trated in Figure !^ 
Moving from the atmospheric grid to the evolutionary grid. Figure 17 of Valenti & Fischer (20051 illustrates that Tor 
the coolest stars under consideration here, which were the focus of their study, variation of metallicity from +0.5 to 
-0.5 dex in [Fe/H] corresponds to a +0.1 to -0.1 dex shift in logg of an evolutionary isochrone. Among hotter stars. 
Figure shows that metallicity uncertainty affects temperatures only minorly, and gravities not at all or minimally. 

We similarly calculated the effect on atmospheric parameter determination when allowing the model color grids to 
vary from +0.5 to -0.5 dex in [M/H], which notably represent the extremes of the m eta llicity range included in our 
sample (less than 1% of stars considered here have |[Fe/H]| > 0.5 dex). Figures 27 & 28 examine in detail the effects 
of metallicity on Teg, log 5 determinations in the relevant uvbyj3 planes. In summary, Jeff variations of up to ^ 1% in 
the {b — y) — ci plane, ^ 2% in the oq — r* plane, and 6 % in the ci — j3 plane are possible with shifts of ± 0.5 dex 
in [M/H]. Notably, however, Tgg variations above the 2% level are only expected in the ci — f3 plane for stars hotter 
than ^ 17000 K, or roughly spectral type B4, of which there are very few in our sample. Similarly metallicity shifts 
of ± 0.5 dex can cause variations of ^ 0.1 dex in logg in the {b — y) — ci and ci — /3 planes, while the same variation 
in the Og — r* plane produces surface gravity shifts closer to ^ 0.05 dex. 


By contrast, metallicity effects are more prominent in color-magnitude techniques. Recently, Nielsen et al. (20131 
executed a Bayesian analysis of the locations in the My vs B — V diagram of Gemini/NICI targets to derive their 
ages including confidence contours for the stellar masses, ages, and metallicities. The work demonstrates correlation 
in this particular color-magnitude diagram of increasing mass and decreasing age with higher metallicity. Metal poor 
stars will have erroneously young ages attributed to them when solar metallicity is assumed. 




Fig. 26.— Comparison of ATLAS9 color grids for different metallicities. 


Confidence Intervals 

All confidence intervals in age and mass quoted in this work are the bounds of the Highest Posterior Density (HPD) 
Region. For a given posterior probability density, p{6\x), the 100(1 — a)% HPD region is defined as the subset, C, of 
9 values: 


c = {e:pie\x)>p*}, 


( 1 ) 


where p* is the largest number such that 


/ 0-p{0\x)>p* 


p{9\x)d9 = 1 — 0 !. 


( 2 ) 


In other words, the HPD region is the set of most probable values (corresponding to the smallest range in 9) that 
encloses 100(1 — o)% of the posterior mass. The HPD method is particularly suited for hnding conhdence intervals of 
skewed probability distributions, such as the stellar age posteriors studied in this work. To find the highest posterior 
density (HPD) region numerically, a function is created that iteratively integrates a normalized posterior PDF above 
a test value of p* while the area/volume under the PDF is less than the desired confidence interval. 

Open Cluster Tables 
Alternative Treatment of Open Clusters 

As described in § |5.2.3| The 1-D marginalized PDF in age for an individual star is computed on a mod el grid that is 


uniformly spaced in log(age). As such, the prior probability of each bin is also encoded in log(age) (see § 5.2.2). Thus 
the resultant PDF is naturally in the units of d p(logr)/dlogr, where p is probability and r is age. 
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Fig . 27.— The effect of metallicity on uvbyjS determinations of temperature, as predicted by model grids of|Castelli Sz Kurucz| (|2006| 
|2004||. In the left-most figure, for given values of Teff, or (b — y)^ the ratio of the temperature given by the grid of metallicity [iVL/tl]=-b.5 to 
tne solar metallicity grid is depicted in the top panel. The bottom panel shows the ratio of the temperature given by a grid of metallicity 
[M/H]=+0.5 to the temperature given by the solar metallicity grid. In the temperature range of interest ( 5 ^ 6500K-8500K, or spectral types 
F5-A4), a shift of 0.5 dex in [M/H] can produce variations up to ~ 1% in Tgff, with the smallest discrepancies occurring at approximately the 
F0-A9 boundary. The middle figure is analogous to the left figure, for the uq — r* grids which are used for stars between 5 ^ 8500K-11000K 
(A3-B9). In this regime, shifts of 0.5 dex in metallicity can produce variations up to ~ 2% in temperature. Finally, for the hottest stars 
(^eff > 11000 K, spectral types B9 and earlier), a 0.5 dex shift in metallicity can produce variations up to ~ 6% in effective temperature. 


{b-y)-c, 





Fig . 28.— The effect of metallicity on uvby^ determinations of surface gravity, as predicted by model grids of|Castelli & Kurucz| (|2006| 
|2004||. In the left-most figure, for given values of logg, or ci, the ratio of the temperature given by the grid of metallicity [M/Jrl]=-U.5 to 
the solar metallicity grid is depicted in the top panel. The bottom panel shows the ratio of the temperature given by a grid of metallicity 
[M/H]=+0.5 to the temperature given by the solar metallicity grid. In the temperature range of interest (ps 6500K-8500K, or spectral 
types F5-A4), a shift of 0.5 dex in [M/H] can produce variations up to ~ 0.1 dex in log^. The middle figure is analogous to the left figure, 
for the ao — r* grids which are used for stars between ^ 8500K-11000K (A3-B9). In this regime, the gravity indicator r* is particularly 
insensitive to metallicity, with shifts of 0.5 dex in metallicity producing variations of only ~ 0.05 dex or less in \ogg. Finally, for the hottest 
stars (Tgff > 11000 K, spectral types B9 and earlier), a 0.5 dex shift in metallicity can produce variations up to ~ 0.1 dex in logg. 

In order to transform p(logr) to p(r) one uses the conversion p(r) = p(logT)/T. Statistical measures other than 
the median^ such as the mean, mode, confidence intervals, etc. will be different depending on whether the PDF being 
quantified is p(logT) or p(r). For example, ^ (r). Strictly speaking, however, both values are meaningful and 

authors frequently choose to report one or the other in the literature. In the case at hand, p(logT) for an individual 
star is more symmetric than the linear counterpart, p(t). As such, one could reasonably argue that 10^'°®'^^ is a more 
meaningful metric than (r). 

In either case, because the PDFs in age or log(age) are both skewed, the median (which, again, is equal regardless 
of whether p(t) or p(logr) is under consideration), is actually the most meaningful quantification of the PDF since it 
is less susceptible to extreme values than either the mean or mode. 

With respect to the open clusters, regardless of whether our analyses are performed in logarthmic or linear space, 
our results favor ages that are younger and older than accepted values for a Per and the Hyades, respectively. 
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TABLE 8 

IC 2602 MEMBERS DEREDDENED uvbyf) PHOTOMETRY AND ATMOSPHERIC PARAMETERS 


Star 

Sp. Type 

ib-y)o 

(mag) 

mo 

(mag) 

CO 

(mag) 

/3 

(mag) 

TeS 

(K) 

logs 

(dex) 

V sin i 
(km s“^) 

HD 91711 

B8 V 

-0.062 

0.146 

0.457 

2.745 

14687 ± 235 

4.467 ± 0.113 

153 

HD 91839 

Al V 

0.025 

0.178 

1.033 

2.904 

9509 ± 152 

4.188 ± 0.091 

146 

HD 91896 

B7 III 

-0.081 

0.093 

0.346 

2.660 

16427 ± 263 

3.782 ± 0.113 

155 

HD 91906 

AO V 

0.016 

0.177 

1.005 

2.889 

9799 ± 157 

4.146 ± 0.113 

149 

HD 92275 

B8 HI/IV 

-0.056 

0.125 

0.562 

2.709 

13775 ± 220 

3.852 ± 0.113 

153 

HD 92467 

B95IH 

-0.026 

0.168 

0.833 

2.851 

11178 ± 179 

4.423 ± 0.113 

no 

HD 92478 

AO V 

0.010 

0.183 

0.978 

2.925 

9586 ± 153 

4.431 ± 0.091 

60 

HD 92535 

A5 V n 

0.104 

0.194 

0.884 

2.838 

8057 ± 129 

4.344 ± 0.145 

140 

HD 92536 

B8 V 

-0.043 

0.131 

0.705 

2.795 

13183 ± 211 

4.423 ± 0.113 

250 

HD 92568 

A M 

0.209 

0.237 

0.625 

2.748 

7113 ± 114 

4.341 ± 0.145 

126 

HD 92664 

B8 III P 

-0.083 

0.118 

0.386 

2.702 

15434 ± 247 

4.145 ± 0.113 

65 

HD 92715 

B9 V nn 

-0.027 

0.136 

0.882 

2.836 

12430 ± 199 

4.362 ± 0.113 

290 

HD 92783 

B85V nn 

-0.033 

0.124 

0.835 

2.804 

12278 ± 196 

4.130 ± 0.113 

230 

HD 92837 

AO IV nn 

-0.007 

0.160 

0.953 

2.873 

10957 ± 175 

4.322 ± 0.113 

220 

HD 92896 

A3 IV 

0.114 

0.193 

0.838 

2.831 

8010 ± 128 

4.425 ± 0.145 

139 

HD 92938 

B3 V n 

-0.075 

0.105 

0.384 

2.690 

15677 ± 251 

4.015 ± 0.113 

120 

HD 92966 

B95V nn 

-0.019 

0.158 

0.930 

2.878 

11372 ± 182 

4.445 ± 0.113 

225 

HD 92989 

A05Va 

0.008 

0.180 

0.982 

2.925 

9979 ± 160 

4.480 ± 0.091 

148 

HD 93098 

Al V s 

0.017 

0.180 

0.993 

2.915 

9688 ± 155 

4.385 ± 0.091 

135 

HD 93194 

B3 V nn 

-0.078 

0.105 

0.357 

2.668 

17455 ± 279 

4.015 ± 0.113 

310 

HD 93424 

A3 Va 

0.060 

0.197 

0.950 

2.890 

8852 ± 142 

4.247 ± 0.113 

95 

HD 93517 

Al V 

0.052 

0.196 

0.976 

2.919 

9613 ± 154 

4.510 ± 0.091 

220 

HD 93540 

B6 V nn 

-0.065 

0.116 

0.476 

2.722 

15753 ± 252 

4.308 ± 0.113 

305 

HD 93549 

B6 V 

-0.066 

0.123 

0.454 

2.729 

15579 ± 249 

4.422 ± 0.113 

265 

HD 93607 

B25V n 

-0.084 

0.102 

0.292 

2.675 

17407 ± 279 

4.098 ± 0.113 

160 

HD 93648 

AO V n 

0.041 

0.188 

1.025 

2.890 

9672 ± 155 

4.157 ± 0.091 

215 

HD 93714 

B2 IV-V n 

-0.092 

0.100 

0.201 

2.647 

18927 ± 303 

3.979 ± 0.113 

40 

HD 93738 

AO V nn 

-0.027 

0.158 

0.842 

2.817 

12970 ± 208 

4.336 ± 0.113 

315 

HD 93874 

A3 IV 

0.071 

0.203 

0.947 

2.896 

8831 ± 141 

4.367 ± 0.091 

142 

HD 94066 

B5 V n 

-0.068 

0.117 

0.439 

2.680 

15096 ± 242 

3.792 ± 0.113 

154 

HD 94174 

AO V 

0.046 

0.193 

0.946 

2.907 

9305 ± 149 

4.391 ± 0.113 

149 
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TABLE 9 

a Persei members dereddened uvbyf ) photometry and atmospheric parameters 


Star 

Sp. Type 

ib-y)o 

(mag) 

mo 

(mag) 

CO 

(mag) 

/3 

(mag) 

Teff 

(K) 

logs 

(dex) 

V sin i 
(km s“^) 

BD+49 868 

F5 V 

0.261 

0.165 

0.459 

2.683 

6693 ± 107 

4.455 

± 

0.145 

20 

HD 19767 

FO V N 

0.176 

0.178 

0.756 

2.765 

7368 ±118 

4.174 

± 

0.145 

140 

HD 19805 

AO Va 

-0.000 

0.161 

0.931 

2.887 

10073 ± 161 

4.344 

± 

0.113 

20 

HD 19893 

B9 V 

-0.031 

0.131 

0.850 

2.807 

12614 ± 202 

4.176 

± 

0.113 

280 

HD 19954 

A9 IV 

0.150 

0.200 

0.794 

2.792 

7632 ± 122 

4.297 

± 

0.145 

85 

HD 20135 

AO P 

-0.011 

0.186 

0.970 

2.848 

10051 ± 161 

3.998 

± 

0.113 

35 

BD+49 889 

F5 V 

0.292 

0.156 

0.418 

2.656 

6430 ± 103 

4.352 

± 

0.145 

65 

BD+49 896 

F4 V 

0.261 

0.168 

0.472 

2.686 

6686 ± 107 

4.410 

± 

0.145 

30 

HD 20365 

B3 V 

-0.079 

0.103 

0.346 

2.681 

16367 ± 262 

4.025 

± 

0.113 

145 

HD 20391 

Al Va n 

0.026 

0.179 

1.006 

2.901 

10415 ± 167 

4.386 

± 

0.091 

260 

HD 20487 

AO V N 

-0.016 

0.151 

0.976 

2.856 

11659 ± 187 

4.198 

± 

0.113 

280 

BD+47 808 

FI IV N 

0.183 

0.179 

0.759 

2.763 

7281 ± 116 

4.062 

± 

0.145 

180 

BD+48 892 

F3 IV-V 

0.246 

0.167 

0.524 

2.696 

6800 ± 109 

4.359 

± 

0.145 

20 

BD+48 894 

FO IV 

0.174 

0.202 

0.734 

2.770 

7416 ± 119 

4.284 

± 

0.145 

75 

HD 20809 

B5 V 

-0.074 

0.109 

0.395 

2.696 

15934 ± 255 

4.097 

± 

0.113 

200 

HD 20842 

AO Va 

-0.005 

0.157 

0.950 

2.886 

10258 ± 164 

4.325 

± 

0.113 

85 

HD 20863 

B9 V 

-0.034 

0.134 

0.810 

2.813 

12154 ± 194 

4.267 

± 

0.113 

200 

BD+49 914 

F5 V 

0.281 

0.170 

0.431 

2.664 

6520 ± 104 

4.395 

± 

0.145 

120 

HD 20919 

A8 V 

0.168 

0.191 

0.757 

2.775 

7463 ±119 

4.259 

± 

0.145 

50 

BD+49 918 

FI V N 

0.186 

0.183 

0.770 

2.755 

7235 ± 116 

3.977 

± 

0.145 

175 

HD 20931 

Al Va 

0.018 

0.174 

0.979 

2.911 

9588 ± 153 

4.342 

± 

0.113 

85 

BD+47 816 

F4 V 

0.271 

0.155 

0.452 

2.672 

6600 ± 106 

4.399 

± 

0.145 

28 

HD 20961 

B95V 

-0.019 

0.163 

0.920 

2.875 

10537 ± 169 

4.344 

± 

0.113 

25 

BD+46 745 

F4 V 

0.274 

0.169 

0.462 

2.674 

6566 ± 105 

4.332 

± 

0.145 

160 

HD 20969 

A8 V 

0.186 

0.192 

0.715 

2.758 

7291 ± 117 

4.239 

± 

0.145 

20 

HD 20986 

A3 V N 

0.046 

0.190 

1.004 

2.896 

9584 ± 153 

4.243 

± 

0.091 

210 

HD 21005 

A5 V N 

0.074 

0.189 

0.987 

2.862 

8266 ± 132 

4.197 

± 

0.145 

250 

HD 21091 

B95IV nn 

-0.019 

0.152 

0.938 

2.856 

12477 ± 200 

4.416 

± 

0.113 

340 

HD 21092 

A5 V 

0.054 

0.218 

0.938 

2.893 

8775 ± 140 

4.311 

± 

0.091 

75 

TYC 3320-1715-1 

F4 V 

0.281 

0.153 

0.469 

2.663 

6495 ± 104 

4.220 

± 

0.145 

no 

HD 21152 

B9 V 

-0.018 

0.158 

0.943 

2.868 

11306 ± 181 

4.353 

± 

0.113 

225 

HD 232793 

F5 V 

0.311 

0.172 

0.377 

2.645 

6274 ± 100 

4.362 

± 

0.145 

93 

HD 21181 

B85V N 

-0.038 

0.122 

0.784 

2.766 

13726 ± 220 

4.119 

± 

0.113 

345 

HD 21239 

A3 V N 

0.045 

0.190 

0.997 

2.910 

9182 ± 147 

4.320 

± 

0.091 

145 

HD 21278 

B5 V 

-0.073 

0.111 

0.398 

2.705 

15274 ± 244 

4.152 

± 

0.113 

75 

HD 21302 

Al V N 

0.022 

0.177 

0.989 

2.888 

10269 ± 164 

4.301 

± 

0.091 

230 

BD+48 923 

F4 V 

0.270 

0.153 

0.464 

2.673 

6603 ± 106 

4.362 

± 

0.145 

20 

HD 21345 

A5 V N 

0.051 

0.208 

0.969 

2.893 

9435 ± 151 

4.324 

± 

0.091 

200 

HD 21398 

B9 V 

-0.030 

0.145 

0.825 

2.837 

11615 ± 186 

4.372 

± 

0.113 

135 

HD 21428 

B3 V 

-0.077 

0.105 

0.363 

2.686 

16421 ± 263 

4.076 

± 

0.113 

200 

HD 21481 

AO V N 

-0.013 

0.164 

0.993 

2.858 

11187 ± 179 

4.141 

± 

0.113 

250 

HD 21527 

A7 IV 

0.093 

0.231 

0.855 

2.856 

8231 ± 132 

4.486 

± 

0.145 

80 

HD 21551 

B8 V 

-0.048 

0.118 

0.673 

2.746 

14869 ± 238 

4.220 

± 

0.113 

380 

HD 21553 

A6 V N 

0.072 

0.206 

0.921 

2.872 

8381 ± 134 

4.414 

± 

0.145 

150 

HD 21619 

A6 V 

0.052 

0.221 

0.935 

2.894 

8843 ± 141 

4.329 

± 

0.091 

90 

BD+49 957 

F3 V 

0.258 

0.168 

0.500 

2.687 

6699 ± 107 

4.334 

± 

0.145 

56 

HD 21641 

B85V 

-0.042 

0.131 

0.721 

2.747 

12914 ± 207 

3.929 

± 

0.113 

215 

BD+49 958 

FI V 

0.198 

0.188 

0.732 

2.739 

7137 ± 114 

3.989 

± 

0.145 

155 

HD 21672 

B8 V 

-0.050 

0.119 

0.649 

2.747 

13473 ± 216 

4.071 

± 

0.113 

225 

BD+48 944 

A5 V 

0.063 

0.220 

0.931 

2.886 

8799 ± 141 

4.305 

± 

0.091 

120 

HD 21931 

B9 V 

-0.029 

0.147 

0.835 

2.829 

11998 ± 192 

4.343 

± 

0.113 

205 
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TABLE 10 

Pleiades members dereddened uvbyf ) photometry and atmospheric parameters 


HD 

Sp. Type 

{b-y)o 

(mag) 

mo 

(mag) 

CO 

(mag) 

/3 

(mag) 

TeS 

(K) 

logg 

(dex) 

V sin i 
(km s“^) 

HD 23157 

A9 V 

0.168 

0.190 

0.725 

2.778 

7463 ± 121 

4.369 

± 

0.145 

100 

HD 23156 

A7 V 

0.111 

0.215 

0.815 

2.837 

8046 ± 130 

4.498 

± 

0.145 

70 

HD 23247 

F3 V 

0.237 

0.174 

0.527 

2.704 

6863 ± 111 

4.424 

± 

0.145 

40 

HD 23246 

A8 V 

0.170 

0.184 

0.758 

2.773 

7409 ± 120 

4.234 

± 

0.145 

200 

HD 23288 

B7 V 

-0.051 

0.120 

0.636 

2.747 

13953 ± 226 

4.151 

± 

0.113 

280 

HD 23302 

B6 III 

-0.054 

0.098 

0.638 

2.690 

13308 ± 216 

3.478 

± 

0.113 

205 

HD 23289 

F3 V 

0.244 

0.164 

0.521 

2.699 

6796 ± 110 

4.387 

± 

0.145 

40 

HD 23326 

F4 V 

0.250 

0.164 

0.514 

2.691 

6741 ± 109 

4.358 

± 

0.145 

40 

HD 23324 

B8 V 

-0.052 

0.116 

0.634 

2.747 

13748 ± 223 

4.126 

± 

0.113 

255 

HD 23338 

B6 IV 

-0.061 

0.104 

0.553 

2.702 

13696 ± 222 

3.772 

± 

0.113 

130 

HD 23351 

F3 V 

0.249 

0.176 

0.507 

2.695 

6755 ± 109 

4.391 

± 

0.145 

80 

HD 23361 

A25Va n 

0.069 

0.201 

0.959 

2.872 

8356 ± 135 

4.309 

± 

0.145 

235 

HD 23375 

A9 V 

0.180 

0.187 

0.710 

2.765 

7336 ± 119 

4.318 

± 

0.145 

75 

HD 23410 

AO Va 

0.004 

0.164 

0.975 

2.899 

10442 ± 169 

4.382 

± 

0.113 

200 

HD 23409 

A3 V 

0.070 

0.202 

0.980 

2.892 

8903 ± 144 

4.270 

± 

0.091 

170 

HD 23432 

B8 V 

-0.039 

0.127 

0.758 

2.793 

12695 ± 206 

4.250 

± 

0.113 

235 

HD 23441 

B9 V N 

-0.029 

0.135 

0.858 

2.822 

11817 ± 191 

4.209 

± 

0.113 

200 

HD 23479 

A9 V 

0.188 

0.166 

0.716 

2.755 

7239 ± 117 

4.212 

± 

0.145 

150 

HD 23489 

A2 V 

0.033 

0.183 

1.012 

2.907 

9170 ± 149 

4.239 

± 

0.091 

no 

HD 23512 

A2 V 

0.057 

0.196 

1.035 

2.909 

8852 ± 143 

4.214 

± 

0.091 

145 

HD 23511 

F5 V 

0.279 

0.174 

0.412 

2.674 

6521 ± 106 

4.477 

± 

0.145 

28 

HD 23514 

F5 V 

0.285 

0.179 

0.443 

2.668 

6450 ± 104 

4.307 

± 

0.145 

40 

HD 23513 

F5 V 

0.278 

0.170 

0.423 

2.673 

6528 ± 106 

4.447 

± 

0.145 

30 

HD 23568 

B95Va n 

-0.024 

0.139 

0.914 

2.847 

11731 ± 190 

4.301 

± 

0.113 

240 

HD 23567 

FO V 

0.159 

0.196 

0.735 

2.788 

7560 ± 122 

4.407 

± 

0.145 

50 

HD 23585 

FO V 

0.168 

0.185 

0.713 

2.780 

7472 ± 121 

4.405 

± 

0.145 

100 

HD 23608 

F5 V 

0.278 

0.177 

0.482 

2.673 

6492 ± 105 

4.185 

± 

0.145 

no 

HD 23607 

FO V 

0.108 

0.203 

0.814 

2.841 

8085 ± 131 

4.534 

± 

0.145 

12 

HD 23629 

AO V 

-0.001 

0.163 

0.986 

2.899 

10340 ± 168 

4.342 

± 

0.113 

170 

HD 23632 

AO Va 

0.006 

0.167 

1.009 

2.899 

10461 ± 169 

4.312 

± 

0.113 

225 

HD 23628 

A4 V 

0.090 

0.189 

0.904 

2.853 

8163 ± 132 

4.381 

± 

0.145 

215 

HD 23643 

A35V 

0.079 

0.194 

0.943 

2.862 

8258 ± 134 

4.301 

± 

0.145 

185 

HD 23733 

A9 V 

0.207 

0.177 

0.672 

2.736 

7066 ± 114 

4.174 

± 

0.145 

180 

HD 23732 

F5 V 

0.258 

0.172 

0.460 

2.688 

6695 ± 108 

4.473 

± 

0.145 

50 

HD 23753 

B8 V N 

-0.046 

0.113 

0.712 

2.736 

13096 ± 212 

3.859 

± 

0.113 

240 

HD 23791 

A9 V-l- 

0.139 

0.214 

0.758 

2.811 

7776 ± 126 

4.480 

± 

0.145 

85 

HD 23850 

B8 III 

-0.048 

0.102 

0.701 

2.695 

13446 ± 218 

3.483 

± 

0.113 

280 

HD 23863 

A8 V 

0.116 

0.201 

0.857 

2.826 

7926 ± 128 

4.354 

± 

0.145 

160 

HD 23872 

Al Va n 

0.032 

0.182 

1.013 

2.894 

10028 ± 162 

4.247 

± 

0.091 

240 

HD 23873 

B95Va 

-0.023 

0.143 

0.907 

2.852 

10897 ± 177 

4.255 

± 

0.113 

90 

HD 23886 

A4 V 

0.068 

0.214 

0.915 

2.880 

8974 ± 145 

4.343 

± 

0.091 

165 

HD 23912 

F3 V 

0.274 

0.154 

0.481 

2.671 

6531 ± 106 

4.242 

± 

0.145 

130 

HD 23924 

A7 V 

0.100 

0.223 

0.852 

2.852 

8121 ± 132 

4.460 

± 

0.145 

100 

HD 23923 

B85V N 

-0.033 

0.124 

0.839 

2.794 

12911 ± 209 

4.159 

± 

0.113 

310 

HD 23948 

Al Va 

0.033 

0.191 

0.984 

2.905 

9237 ± 150 

4.307 

± 

0.091 

120 

HD 24076 

A2 V 

0.008 

0.168 

0.923 

2.867 

10196 ± 165 

4.298 

± 

0.091 

155 

HD 24132 

F2 V 

0.245 

0.149 

0.597 

2.692 

6744 ± 109 

4.182 

± 

0.145 

230 
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TABLE 11 

Hyades members dereddened uvbyjS photometry and atmospheric parameters 



HD 

Sp. Type 

{b-y)o 

mo 

CO 

/3 

T'eff 


logff 

V sin i 




(mag) 

(mag) 

(mag) 

(mag) 

(K) 


(dex) 

(km s“^) 

HD 

26015 

F3 V 

0.252 

0.174 

0.537 

2.693 

6732 

± 

109 

4.244 

± 

0.145 

25 

HD 

26462 

FI IV-V 

0.230 

0.165 

0.596 

2.710 

6916 

± 

112 

4.291 

± 

0.145 

30 

HD 

26737 

F5 V 

0.274 

0.168 

0.477 

2.674 

6558 

± 

106 

4.263 

± 

0.145 

60 

HD 

26911 

F3 V 

0.258 

0.176 

0.525 

2.690 

6682 

± 

108 

4.228 

± 

0.145 

30 

HD 

27176 

A7 m 

0.172 

0.187 

0.785 

2.767 

7380 

± 

120 

4.087 

± 

0.145 

125 

HD 

27397 

FO IV 

0.171 

0.194 

0.770 

2.766 

7410 

± 

120 

4.173 

± 

0.145 

100 

HD 

27429 

F2 VN 

0.240 

0.171 

0.588 

2.693 

6828 

± 

111 

4.270 

± 

0.145 

150 

HD 

27459 

FO IV 

0.129 

0.204 

0.871 

2.812 

7782 

± 

126 

4.198 

± 

0.145 

35 

HD 

27524 

F5 V 

0.285 

0.161 

0.461 

2.656 

6461 

± 

105 

4.213 

± 

0.145 

110 

HD 

27561 

F4 V 

0.270 

0.162 

0.482 

2.677 

6594 

± 

107 

4.284 

± 

0.145 

30 

HD 

27628 

A2 M 

0.133 

0.225 

0.707 

2.756 

7944 

± 

129 

4.743 

± 

0.145 

30 

HD 

27819 

A7 IV 

0.080 

0.209 

0.982 

2.857 

8203 

± 

133 

4.170 

± 

0.145 

35 

HD 

27901 

F4 V N 

0.238 

0.178 

0.597 

2.704 

6837 

± 

111 

4.233 

± 

0.145 

110 

HD 

27934 

A5 IV-V 

0.064 

0.201 

1.053 

2.867 

8506 

± 

138 

3.884 

± 

0.091 

90 

HD 

27946 

A7 V 

0.149 

0.192 

0.840 

2.783 

7584 

± 

123 

4.112 

± 

0.145 

210 

HD 

27962 

A3 V 

0.020 

0.193 

1.046 

2.889 

9123 

± 

148 

4.004 

± 

0.091 

30 

HD 

28024 

A9 IV- N 

0.165 

0.175 

0.947 

2.753 

7279 

± 

118 

3.503 

± 

0.145 

215 

HD 

28226 

A M 

0.164 

0.213 

0.771 

2.775 

7493 

± 

121 

4.248 

± 

0.145 

130 

HD 

28294 

FO IV 

0.198 

0.173 

0.694 

2.745 

7174 

± 

116 

4.194 

± 

0.145 

135 

HD 

28319 

A7 III 

0.097 

0.198 

1.011 

2.831 

7945 

± 

129 

3.930 

± 

0.145 

130 

HD 

28355 

A7 m 

0.112 

0.226 

0.908 

2.832 

7930 

± 

128 

4.207 

± 

0.145 

140 

HD 

28485 

FO V-l- N 

0.200 

0.192 

0.717 

2.740 

7129 

± 

115 

4.035 

± 

0.145 

150 

HD 

28527 

A5 m 

0.085 

0.218 

0.964 

2.856 

8180 

± 

133 

4.194 

± 

0.145 

100 

HD 

28546 

A7 m 

0.142 

0.234 

0.796 

2.809 

7726 

± 

125 

4.354 

± 

0.145 

30 

HD 

28556 

FO IV 

0.147 

0.202 

0.814 

2.795 

7645 

± 

124 

4.244 

± 

0.145 

140 

HD 

28568 

F5 V 

0.274 

0.168 

0.466 

2.676 

6564 

± 

106 

4.315 

± 

0.145 

55 

HD 

28677 

F2 V 

0.214 

0.176 

0.654 

2.725 

7032 

± 

114 

4.161 

± 

0.145 

100 

HD 

28911 

F5 V 

0.283 

0.163 

0.459 

2.663 

6481 

± 

105 

4.249 

± 

0.145 

40 

HD 

28910 

A9 V 

0.144 

0.200 

0.830 

2.796 

7659 

± 

124 

4.213 

± 

0.145 

95 

HD 

29169 

F2 V 

0.236 

0.183 

0.567 

2.708 

6880 

± 

111 

4.321 

± 

0.145 

80 

HD 

29225 

F5 V 

0.276 

0.171 

0.461 

2.675 

6547 

± 

106 

4.316 

± 

0.145 

45 

HD 

29375 

FO IV-V 

0.187 

0.187 

0.740 

2.754 

7257 

± 

118 

4.106 

± 

0.145 

155 

HD 

29388 

A5 IV-V 

0.062 

0.199 

1.047 

2.870 

8645 

± 

140 

3.927 

± 

0.091 

115 

HD 

29499 

A M 

0.140 

0.231 

0.826 

2.810 

7713 

± 

125 

4.266 

± 

0.145 

70 

HD 

29488 

A5 IV-V 

0.080 

0.196 

1.017 

2.852 

8127 

± 

132 

4.025 

± 

0.145 

160 

HD 

30034 

A9 IV- 

0.150 

0.195 

0.813 

2.791 

7610 

± 

123 

4.218 

± 

0.145 

75 

HD 

30210 

A5 m 

0.091 

0.252 

0.955 

2.845 

8126 

± 

132 

4.181 

± 

0.145 

30 

HD 

30780 

A9 V-l- 

0.122 

0.207 

0.900 

2.813 

7823 

± 

127 

4.141 

± 

0.145 

155 

HD 

31845 

F5 V 

0.294 

0.165 

0.439 

2.658 

6396 

± 

104 

4.229 

± 

0.145 

25 

HD 

32301 

A7 IV 

0.079 

0.202 

1.034 

2.847 

8116 

± 

131 

3.975 

± 

0.145 

115 

HD 

33254 

A7 m 

0.132 

0.251 

0.835 

2.824 

7797 

± 

126 

4.306 

± 

0.145 

30 

HD 

33204 

A7 m 

0.149 

0.245 

0.803 

2.796 

7634 

± 

124 

4.270 

± 

0.145 

30 

HD 

25202 

F4 V 

0.206 

0.172 

0.695 

2.724 

7082 

± 

115 

4.064 

± 

0.145 

160 

HD 

28052 

FO IV-V N 

0.153 

0.183 

0.934 

2.767 

7431 

± 

120 

3.733 

± 

0.145 

170 

HD 

18404 

F5 IV 

0.269 

0.169 

0.481 

2.680 

6605 

± 

107 

4.299 

± 

0.145 

0 

HD 

25570 

F4 V 

0.249 

0.147 

0.557 

2.688 

6752 

± 

109 

4.183 

± 

0.145 

34 

HD 

40932 

A2 M 

0.079 

0.205 

0.978 

2.853 

8224 

± 

133 

4.191 

± 

0.145 

18 


TABLE 12 

Open Cluster Ages: Analysis in Linear Age 

Cluster 

Lit. Age 
(Myr) 

Models 


Summed PDF 
Median 
(Myr) 

Summed PDF 
68% C.I. 
(Myr) 

Multiplied PDF 
Median 
(Myr) 

Multiplied PDF 

68% C.I. 

(Myr) 

IC 2602 

46l| 

Ekstrom et al. 

(2012 

1 22 

3-39 

41 

41-42 



Bressan et al.l 

2012 

24 

3-40 

40 

37-43 

ol Persei 

Ekstrom et al. 

(2012 

I 41 

3-68 

63 

61-68 


125^« 

Bressan et ai.l 

2012 

45 

3-71 

62 

58-66 

Pleiades 

Ekstrom et al. 

(2012 

1 61 

3-113 

125 

122-131 



Bressan et al.l 

2012 

77 

3-117 

112 

107-120 

Hyades 

Ekstrom et al. 

(2012 

I 118 

3-403 

677 

671-690 



Bressan et al.| 


288 

17-593 

738 

719-765 


Literature ages (column 2) come from the sources referenced in §[^ For each set of evolutionary models, the median and 68% confidence 
interval are computed for both the summed PDF (columns 4,5) and multiplied PDF (columns 6,7). Note, the Hyades analysis includes 
the blue straggler HD 27962 and the spectroscopic binary HD 27268. Excluding these outliers results in a median and 68% confidence 
interval of 322 Myr [17-650 Myr] of the summed PDF or 784 Myr [749-802 Myr[ of the multiplied PDF, using the B12 models. 
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Fig. 29.— Left panels: ID marginalized, normalized posterior PDFs in age, calculated from|Bressan et al.|||2012|| evolutionary models, for 
individual open cluster members. Black, teal, and red histograms represent early, intermediate, and late group stars, respectively. Middle 
panels: Sums of the individual PDFs depicted on the left. This figure shows the total probability associated with the 200 age bins between 
log(age/yr)=6.5 to 10. The grey shaded regions indicate the currently accepted ages of IC 2602 (46^5 Myr), a Per (90dil0 Myr), the 
Pleiades (125di8 Myr), and the Hyades (625di50 Myr). Right panels: Products of the individual PDFs depicted in the left panels. The 
grey shaded regions again depict the accepted literature age ranges of each cluster. 
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